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ABSTRACT 

This  document  outlines  the  Helicopter  Slung  Load  Simulation  (HSLSIM)  Toolbox  which  is  a  set  of 
utilities  for  the  simulation,  analysis  and  display  of  the  flight-dynamic  response  of  helicopters  with 
various  external  load  configurations  within  the  MATLAB®  software  environment.  Instructions 
and  examples  are  provided  for  its  operation  and  subsequent  modifications.  The  helicopter 
studied  is  the  CH-47D,  with  the  load  types  including  rectangular  and  cylindrical  containers  as 
well  as  plate  and  airfoil  shapes. 
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Helicopter  Slung-Load  Simulation  Toolbox 
for  use  with  MATLAB® 


Executive  Summary 

The  motivation  behind  this  task  comes  from  a  requirement  by  the  Austrahan  Army  to 
estabhsh  the  safe  operating  Limits  of  hehcopters  with  externally  slimg  loads  in  flight.  The 
dynamic  behaviour  of  slimg-loads  can  be  extremely  difficult  to  predict  and  depends  on 
such  things  as  the  load  shape  and  density^  the  hehcopter  airspeed  and  control  inputs  as 
weU  as  the  shng  configuration.  Many  incidents  have  been  reported  where  the  load  has 
oscillated  uncontrollably  or  completely  diverged  due  to  aerodynamic  excitation  resulting 
in  damage  and/or  loss  of  equipment.  The  interaction  of  these  elements  is  not  well 
understood.  Hence,  a  major  requirement  of  the  study  was  to  develop  a  comprehensive 
hehcopter  slimg-load  model  to  provide  a  better  imderstanding  of  the  system  dynamics 
and  various  effects  involved.  Since  the  initial  development  the  model  has  been  refined, 
extended  and  packaged  together  in  the  form  of  a  Toolbox  that  is  accessed  from  within  the 
MATLAB®  numerical  computing  software  environment. 

This  document  outlines  the  Hehcopter  Slung-Load  Simulation  (HSLSIM)  software  package 
and  provides  instructions  and  examples  for  its  operation  and  subsequent  modification. 
The  graphical  user  interface  is  presented  along  with  the  various  configuration  dialogs  and 
output  windows  generated  during  each  stage  of  typical  operation.  These  dialogs  include 
configuration  options  for  the  hehcopter,  loads,  control  inputs  and  initial  conditions.  The 
three  main  stages  of  operation  are:  Balance,  in  which  the  load(s)  are  balanced  in  the  free- 
stream  below  the  fixed  hehcopter;  Trim,  where  the  hehcopter  and  loads  are  both  trimmed 
at  the  flight  state  specified  and;  Simulation,  in  which  the  hehcopter-load  system  response 
to  specified  control  inputs  is  simulated  over  the  time  range.  The  simulation  may  then  be 
replayed  either  from  within  MATLAB®  or  in  an  external  Virtual  Reahty  Modelling 
Language  (VRML)  browser.  As  an  example,  the  trim  and  simulation  procedure  for  a  CH- 
47D  hehcopter  carrying  a  3000L  water  tank  is  given.  Using  the  graphical  user  interface,  the 
hehcopter-load  system  is  trimmed  and  flown  at  a  nominal  speed  (chosen  as  100  ft/s)  with 
a  lateral  doublet  input.  Then,  using  a  simple  script  to  automate  the  process,  the  same 
configuration  is  flown  over  a  range  of  airspeeds  and  load  masses.  Two  methods  for 
presentation  of  the  resulting  data  are  then  shown. 

With  the  Hehcopter  Slimg-Load  Simulation  Toolbox,  the  Austrahan  Army  whl  be  able  to 
gain  further  insight  into  the  dynamic  behaviour  of  various  loads  under  different  flight 
conditions.  The  user  can  examine  the  characteristics  of  particular  hehcopter-load 
configurations,  or  assess  the  performance  of  such  configurations  over  a  range  of  different 
flight  states. 
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1.  Introduction 

The  Helicopter  Sltmg-Load  Simiilation  package  (HSLSIM)  has  been  developed  at  the  Defence 
Science  and  Technology  Organisation  (DSTO)  as  a  fundamental  part  of  Task  ARM  07/038  and 
preceding  tasks  sponsored  by  Australian  Army  Aviation.  In  the  pash  the  operations  of 
helicopters  carrying  externally  slung  loads  has  often  been  limited  and^  in  some  cases^ 
seriously  hindered  by  stabihty  and  control  problems.  A  program  was  consequently  initiated 
within  DSTO  to  use  computer  modelling  and  simulation  to  assist  in  defining  the  operational 
hmits  of  the  Australian  Army  Chinook  CH-47D  helicopter  when  carrying  slung  loads.  The 
first  phase  in  this  program  entailed  the  development  of  a  helicopter  slimg-load  model  for 
simulation  and  analysis  in  order  to  provide  a  better  imderstanding  of  the  system  dynamics 
and  various  effects  involved. 

HSLSIM  was  written  in  M ATL AB'®^  a  high-performance  numerical  computing  environment 
that  is  built  arotmd  matrix  mathematics^  and  therefore  amenable  to  dynamic  modelling  and 
simulation  work.  MALTA B®  also  provides  an  ideal  workplace  in  which  various  analyses  can 
be  conducted  following  the  simulation  itself. 

2.  System  Description 

The  following  section  provides  a  general  description  of  the  system^  including  the  model 
formulation  and  implementation.  For  a  more  detailed  description^  refer  to  the  previous  report 
by  Stuckey  [1]. 

2.1  Model  Formulation 

Hehcopter  slimg-load  systems  fall  into  a  class  of  multibody  systems  approximated  by  two  or 
more  rigid  bodies  connected  by  massless  hnks.  The  huks  can  be  considered  either  elastic  or 
inelastic^  although  the  rigid-body  assumption  excludes  any  helicopter  or  load  elastic  modes. 
Typically^  the  system  is  characterised  by  the  configuration  geometry^  mass^  inertia^  and 
aerodynamic  behaviour  of  both  hehcopter  and  load,  as  weU  as  the  elastic  properties  of  the 
hnks. 

In  general  terms,  the  system  of  interest  consists  of  a  single  hehcopter  supporting  one  or  more 
loads  by  means  of  some  suspension.  The  model  is  comprised  of  n  rigid  bodies,  with  m 
straight-hne  hnks  supporting  a  single  force  in  the  direction  of  the  hnk.  If  the  hnks  are 
modehed  as  inelastic,  c<m  constraints  are  imposed  on  the  motion  of  the  bodies  and  the 
system  has  d  ^  n^'  6  -  c  degrees-of-freedom  (dof).  h  the  hnks  are  modehed  as  elastic,  there  are 
;7*6  dof. 

A  number  of  simplifying  assumptions  have  been  made  in  the  model.  These  include  the 
exclusion  of  cable  mass,  cable  aerodynamics  and  r o tor-do wnw ash  effects.  Despite  these 
limitations,  the  system  has  proven  adequate  for  simulation  studies  [2]  in  which  the  low- 
frequency  behaviour  is  of  primary  interest  and  the  hehcopter  is  initiahy  trimmed  in  forward 
flight. 
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The  simidation  model  used  is  based  on  the  helicopter  slimg-load  system  introduced  by 
Cicolani^  et  al  [3].  In  this  formulation^  the  general  system  equations  of  motion  are  obtained 
from  the  Newton- Euler  equations  in  terms  of  generalised  coordinates  and  velocities. 
Following  the  expbdt  constraint  method^  which  utilises  d'Alembert's  principle^  the  system  is 
partitioned  into  coordinates  such  that  the  motion  due  to  cable  stretching  is  separated  from 
that  due  to  rigid-body,  coupled  dynamics.  As  a  consequence,  the  constraint  forces  appear 
explicitly  and  a  solution  to  the  resultant  generalised  accelerations  is  determined  by  assuming 
a  simple  spring  model  for  the  cable. 

It  is  also  possible  to  obtain  a  solution  to  the  inelastic  approximation  by  nulling  the  stretching 
coordinates  to  obtain  an  expUcit  relation  for  the  suspension  forces.  The  result  is 
computationally  more  efficient  than  conventional  procedures  and  is  readily  integrated  with 
the  formulation  for  elastic  suspension.  Another  benefit  of  the  formulation  is  that  it  is  easily 
applied  to  complex,  multiple  body  systems,  as  in  the  software  package  presented.  Aside  from 
the  core  heUcopter  model,  all  code  development  has  been  done  in  the  MATLAB®  [4]  numerical 
computing  environment,  which  provides  a  high-performance  language,  amenable  to 
modelling  and  simulation  type  work. 

2.2  Implementation  in  MATLAB® 

The  Helicopter  Slimg-Load  Simulation  program  HSLSIM  consists  of  several  modules,  written 
in  the  MATLAB®  language.  These  include  the  main  script,  an  optimisation  routine,  a 
differential  equation  solution,  an  integration  function,  several  flight-dynamic  models,  and 
various  output  and  replay  functions.  There  is  also  a  graphical  user  interface  for  simplified 
control  of  the  primary  program  functions.  Alternatively,  the  simulation  can  be  run  through  a 
main  script,  which  generates  the  control  inputs,  configures  the  heHcopter-load  system 
properties  (geometric  and  inertial),  sets  the  initial  system  state,  and  then  executes  the  trim  and 
integration  functions.  The  integration  function  is  problem  independent  and  based  on  an 
algorithm  which  combines  various  order  Runge-Kutta  formulas  for  the  solution  of  ordinary 
differential  equations.  It  requires  an  integration  function  tailored  to  the  problem  at  hand, 
which  provides  a  point  solution  to  the  differential  equation.  For  the  helicopter  slung-load 
simulation,  this  function  represents  the  core  of  the  code  and  implements  the  aerodynamic 
models  for  both  helicopter  and  loads. 

Initially,  an  approximated  model  of  the  CH-47B,  named  HSL,  and  based  on  a  set  of  linearised 
state-space  representations  at  several  airspeeds  was  used.  This  model  was  derived  from  a  full 
nonlinear  simulation  model  developed  by  the  Boeing  Vertol  Company  and  later  adapted  for 
use  at  the  NASA  Ames  Research  Center  [5]  [6] .  However,  the  model  was  inadequate  for  large 
manoeuvres  or  long  simulations,  so  a  higher  fidelity  model  was  introduced.  This  model, 
named  ROTORGEN,  was  developed  by  Heffley  [7]  for  the  US  Army  Aeroflightdynamics 
Electorate  imder  a  NASA  contract  to  Hoh  Aeronautics  Inc.  It  is  described  as  a  "minimal- 
complexity  generic  rotorcraft  model"  intended  for  maimed  simulation  of  large  military 
helicopters  and,  in  particular,  the  CH-47D  Chinook  tandem  rotor  helicopter,  though  it  is  also 
capable  of  simulating  CH-53D  and  UH-IH  helicopters.  ROTORGEN  implements  a  rotor 
inflow  model  based  on  Glaurt's  representation  of  thrust,  with  the  orientation  (incidence)  of  the 
tip  path  plane  defined  by  a  set  of  flapping  equations  and  the  body  forces  based  on  a  quadratic 
fluid-dynamics  formulation,  applicable  to  low-speed  flight.  The  model  is  actually  a 
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combination  of  two  existing  flight  models:  the  Extended  Stability  Derivative  (ESD)  model 
developed  for  NASA^  and  the  ROTORGEN  thrust  model  developed  for  the  US  Army.  As 
such^  the  ROTORGEN  model  has  a  modular  structure^  which  combines  several  features  of  the 
original  ESD  model.  These  include  a  primary  Fhght  Control  System  (FCS)^  rotor  and  body 
forces^  groimd  effects^  simple  sltmg  load  dynamics^  and  a  Stabhity  and  Control  Augmentation 
System  (SCAS).  HSL  was  written  in  M  ATLAB®,  and  easily  integrated  into  the  simulation  code 
structure.  ROTORGEN^  on  the  other  hand^  was  written  in  Fortran- 77  and  therefore  required 
an  external  (MEX)  routine  to  interface  with  the  MATLAB'®  environment  and  the  simulation 
code.  It  appears  as  a  compiled^  Dynamic  Link  Library  (DLL)  file  in  the  base  directory. 

HSLSIM  also  offers  the  choice  of  several  generic  loads  including  a  box^  cylinder^  plate  and 
airfoil  shapes.  The  box  shape  was  derived  from  a  MIL  VAN  container  —  a  common  hehcopter 
cargo  used  in  many  commercial  and  mhitary  operations.  The  dimensions  of  a  MIL  VAN 
container  are  20  x  8  x  8  ft  and  the  mass  typicaUy  varies  from  4000  lb  (empty)  to  20000  lb  (fuh). 
The  aerodynamic  model  for  the  MILVAN  was  developed  by  Ronen  [8]  and  constitutes  a 
combination  of  aerodynamic  models  obtained  from  several  pubhshed  wind-tunnel  test 
results^  including  unsteady  aerodynamic  effects.  However^  the  actual  model  used^  whhe 
comprehensive^  is  based  purely  on  static  data.  The  cylinder  and  plate  shapes  are  all  compiled 
from  ESDU  data  sheets  [9]  and  Hoerner  [10]  [11]^  while  the  airfoil  -  a  finite-span  NACA  0012 
profile  -  is  extracted  from  the  Comprehensive  Analytical  Model  of  Rotor  craft  Aerodynamics 
and  Dynamics  (CAMRAD/ JA)  code  [12]. 

Normally^  the  simulation  model  is  first  initialLsed  with  a  set  of  flight  conditions  and  given  a 
pre-determined  sequence  of  control  inputs  to  be  used  in  the  computation  of  a  response.  The 
model  is  then  trimmed  according  to  the  specified  flight  conditions  imtil  an  equifibrium  state  is 
reached,  and  then  integrated  over  the  time  range  to  produce  a  set  of  simulation  data.  The  data 
comprises  the  hehcopter  and  load  positions,  orientations,  the  control  inputs  and  geometric 
information  regarding  the  hehcopter-load  configuration.  As  for  most  aeronautical  systems 
[13],  the  positions  are  specified  in  right-handed  Cartesian  axes  (with  positive-z  downward) 
and  the  orientations  in  Euler  angles. 

In  addition  to  the  trim  and  simulation  software,  a  fimction  for  the  replay  of  simulations  was 
written.  This  function,  REPLAY,  displays  rudimentary  models  of  the  hehcopter  and  load(s)  in 
3  dimensions,  with  the  controls  in  a  separate  window.  A  menu  window  ahows  the  user  to 
start,  stop,  pause  and  resume  the  simulation  replay,  as  w eh  as  providing  several  options  for 
the  display. 

The  simulation  output  can  also  be  exported  to  a  file  in  the  Virtual  Reahty  Modelling  Language 
(VRML)  format.  Using  a  VRML  browser,  or  a  WWW  browser  plug-in,  such  as  the  Cortona 
VRML  plug-in  [14],  the  simulation  can  then  be  replayed  with  higher  detail  in  a  similar  3- 
dimensional  space  in  real  time.  Furthermore,  the  file  can  be  compressed  into  a  compact  binary 
format  using  GZIP  [15],  which  can  stih  be  interpreted  by  the  VRML  browser,  but  is  smah 
enough  for  fast  transmission  over  the  Internet. 
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3.  User's  Guide 


3.1  Getting  Started 

The  graphical  user  interface  GHSLSIM  provides  simplified  configuration  and  control  of  the 
heHcopter-load  system  for  trim  and  simulation.  Many  options  are  automated  and  limited  from 
the  full  capabihty  of  the  program^  which  is  discussed  in  Section  3.7. 

The  program  provides  various  configuration  options,  as  well  as  controls  for  balance,  trim  and 
simulation  stages.  There  are  also  controls  for  the  output  and  display  of  the  simulation,  once 
rim.  Balance  refers  to  the  determination  of  an  equilibrium  state  for  the  load(s)  only,  slimg 
beneath  a  fixed  helicopter.  This  intermediate  step  is  performed  to  attain  a  state  closer  to  the 
actual  trim  state  than  the  default,  initial  one  specified  and  is  usually  conducted  prior  to  the 
trim  stage.  Since  the  helicopter  is  fixed  during  this  stage,  it  will  generally  not  be  sufficiently 
trimmed  following  balance.  Therefore,  it  is  necessary  to  trim  the  helicopter  andload(s)  using 
the  trim  control,  regardless  of  whether  it  has  been  balanced  or  not. 

GHSLSIM  can  be  started  by  executing  the  following  command  from  inside  the  MATLAB® 
workspace: 

»  ghslslm 

If  a  configuration  file  exists,  it  may  be  loaded  automatically  by  adding  its  name  after  the 
command.  For  example,  executing: 

»  ghslslm  ghs1deino_3k1t 

will  load  the  configuration  options  stored  within  gh  sl  demo_3 kl  t .  mat.  If  no  file  is  specified, 
the  default  configuration  file,  ghs1  default. mat,  is  used. 

3.2  Main  Interface 

The  main  interface  is  shown  in  Figure  3.1.  There  is  a  list  of  commands  along  the  top  menu  bar 
as  well  as  a  set  of  fields  and  buttons  that  constitute  the  simulation  configuration  and  control 
interface.  The  File  menu  commands  are  accessed  through  the  top  menu  only,  but  the 
Configuration,  Sinnulation  and  Output  commands  can  all  be  accessed  either  through  the  top  menu, 
or  on  the  main  interface  panel. 
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Figure  3.1\  Simulation  Configuration  Interface 


3.2,1  File  Menu  Commands 

As  well  as  specifying  a  configuration  file  on  the  command  line^  it  is  possible  to  load  a 


configuration  file  from  the  file  menu  via  the  command:  |File|  Open.  This  loads  the 


configuration  details  and  the  simulation  output  data,  if  present,  into  the  program  workspace. 
If  a  file  with  simulation  data  is  loaded,  then  the  user  may  execute  any  of  the  Output 
commands  on  the  main  interface  without  having  to  re-run  the  simulation. 

In  order  to  save  a  configuration,  or  both  the  configuration  and  any  simulation  output  data,  the 
user  simply  creates  a  new  file,  or  selects  an  existing  one  from  the  File  -♦  Save  or  the 


File 


Save  As  menu  commands.  If  no  file  has  been  specified,  the  default  file. 


ghslde fault , mat  in  the  current  working  directory  is  used. 


With  the  Qpen|  command,  the  program  expects  a  full  configuration  file  with  all  components 
present,  and  will  produce  a  warning  otherwise.  Similarly,  saving  to  file  using  the  Save  or 
Save  As  commands  saves  all  of  the  configuration  information  and  simulation  data,  if  present. 
If  the  user  wishes  to  load  or  save  only  one  component  of  the  configuration,  then  they  must  use 


the  import  or  export  commands,  respectively.  The  |File|  Import  Config  Data  command  will 


import  whatever  components  are  present  in  the  chosen  partial-configuration  file.  All  other 


components  will  remain  as  they  were,  prior  to  importing.  The  File!  [Export  Config  Data|  has 


several  sub-commands:  Simulation,  Helicopter,  Load,  Control,  and  Trim.  These  can  be  used  to 
export  single  components  from  the  configuration,  which  can  later  be  imported  as  detailed 
above. 
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Preferences. 


These  include 


Preferences  for  the  pro  gram  are  set  in  the  main  menu  via:  [pil^  ■ 
three  fields  for  the  location  of  a  VRML-capable  VJeh  Browser,  the  VP^I'AL  Directory  and  a 
GZIP  Executable,  as  shown  in  Figure  3.2. 


Figure  3,2\  Program  Preferences  Interface 

If  not  present  or  incorrect  the  Web  Browser  location  can  be  typed  directly  into  the  text  fielct  or 
located  by  simply  clicking  the  button  and  manually  searching  through  the  Windows  directoiy 
structure.  The  VRML  Directory  is  the  location  of  all  heUcopter  and  load  VRML  models 
included  in  the  distribution,  and  the  directory  into  which  any  generated  VRML  simulations 
will  be  written.  This  string  ends  with  a  trailing  directory  delimiter.  If  the  GZJP  Executable  is 
installed,  the  last  field  may  be  filled,  or  located  as  with  the  Web  Browser. 

The  preferences  also  include  four  checkboxes  for  various  windows  and  screen  output  that  can 
be  activated  during  balance,  trim,  and  simulation.  These  are  the  Replay  Window,  the  Variable 
Tracking  Windorw,  the  Counter  Window,  and  the  Verbose  Output  flag  The  first  two  will  be 
explained  in  more  detail  in  Section  3.4.  The  Counter  Window  is  just  a  visual  indicator  of  the 
progress  during  simulation  and  post-processing  stages.  Checking  the  Verbose  Output 
preference  will  cause  the  program  to  display  numerical  data  during  trim  and  simulation 
iteration  loops. 

3.2.2  Simulation  Configuration 

Each  of  the  fields  and  options  displayed  in  the  main  interface  are  required  by  the  simulation. 

The  Simulation  time  tN,  Integration  timestep  dt,  and  Number  of  steps  N  are  aU  related  via  the 
following  equation: 


tj^=dt^N  (1) 

Altering  the  Simulation  time  will  also  adjust  the  Number  of  steps  and  altering  the  Integration 
timestep  or  the  Number  of  steps  will  adjust  the  Simulation  time,  such  that  aU  parameters  conform 
to  the  above  relationship.  In  general,  a  smaller  timestep  wUl  yield  more  accurate  simulation 
results  and  reduce  the  possibUity  of  numerical  error  and  subsequent  divergence.  For  most 
configurations,  a  value  between  0.05  s  and  0.1  s  (20  Hz  and  10  Hz)  should  be  sufficient 
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The  next  field  defines  the  Number  ofloads^  which  can  be  any  integer  between  zero  and  three; 
zero  referring  to  the  helicopter-only  configuration.  The  helicopter  attachment  point  (hook)  to 
which  each  load^s  slings  are  attached  is  automatically  determined.  For  imusual  sling 
configurations^  the  user  must  write  a  script  as  detailed  in  Section  3.7.2^  to  construct  fhe  model 
and  rim  the  trim  and  simulation. 

The  following  three  fields  set  the  Trim  airspeed^  in  each  inertial  axis^  at  which  the  helicopter 
will  be  trimmed.  These  are  (in  ft/s)  the  Forward  speedy  rate  of  Descent^  and  Sideslip,  It  is  not 
possible  to  trim  the  aircraft  with  a  constant  (non- zero)  angular  rate^  such  as  in  a  steady  turn. 
The  last  field  sets  the  trim  altitude  of  the  helicopter^  in  ft. 

Options  presented  in  the  interface  comprise  the  Simulation  axes,  Angular  representation, 
Simulation  type.  Sling  model.  Trim  routine  and  Simulation  routine.  The  Simulation  axes  define 
which  axes  are  to  be  used  in  the  trim  and  simulation.  Both  longitudinal  and  lateral  axes 
constitute  a  reduced  subspace  in  which  the  system  is  free  to  respond  (reduced  degrees-of- 
freedom).  The  longitudinal  axes  include  surge^  heave  and  pitch  motion^  while  the  lateral  axes 
include  heave^  sway^  roll  and  yaw  motion.  The  combined  axes  represent  a  fuU^  imconstrained 
system  with  the  number  of  degrees-of-freedom  equal  to  n^6  for  an  elastic  sling  configuration^ 
where  n  denotes  the  number  of  bodies  (helicopter  +  loads).  Unless  there  is  a  specific  need  to 
excite  only  fhe  longitudinal  or  lateral  modes^  it  is  recommended  that  the  combined  axes  be 
used. 

The  angular  displacements  and  rates  for  the  system  can  be  represented  using  either  Euler 
angles  or  quaternions.  However^  this  only  applies  to  the  simulation  stage^  as  the  helicopter  and 
loads  can  only  be  balanced  and  trimmed  using  Euler  angles.  For  most  manoeuvres^  an  Euler 
representation  is  adequate. 

If  the  Control  System  Toolbox  [16]  is  present  within  MATLAB®^  it  is  possible  to  simulate  the 
system  response  using  a  linear  approximation^  which  is  much  faster  than  the  nonlinear 
solution.  These  constitute  the  two  Simulation  types  available  to  the  user.  However^  since  the 
linearised  model  is  obtained  by  computing  a  Jacobian  at  trim^  it  is  inaccurate  for  anything 
other  than  very  short  runs.  In  facb  given  the  high  nonlinearity  of  the  coupled  helicopter-load 
system^  a  linear  simulation  is  reaUy  only  suitable  for  helicopter-only  configurations.  If  the 
Control  System  Toolbox  is  not  present^  a  nonlinear  simulation  must  be  rvm. 

The  two  Sling  models  that  can  be  employed  are  inelastic  and  elastic,  each  incorporating  a 
significantly  different  solution  to  the  multi-body  equations  within  the  balance^  trim  and 
simulation  stages.  The  inelastic  model  uses  an  explicitly  constrained  representation  and 
therefore^  has  less  dof  than  its  elastic  equivalent.  Typically^  the  variables  used  in  an  inelastic 
model  include  the  helicopter  position  and  orientation^  the  orientation  of  one  or  more  cables 
and  the  orientation  of  the  load(s).  The  variables  used  in  an  elastic  model  include  the  position 
and  orientation  of  both  helicopter  and  load(s).  Generally  speakings  less  dof  means  a  lower 
complexity  and  since  the  only  matrix  inversion  required  for  solution  is  calculated  explicitly,  it 
also  means  a  faster  solution.  In  addition,  the  inelastic  model  has  a  higher  degree  of 
orthogonality,  and  is  consequently  much  better  suited  for  trim  optimisation.  This  will  be 
explored  in  the  examples  that  foUow  (Sections  3.6  &  3.7).  Note  that  during  balance  and  trim,  it 
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is  possible  to  switch  between  these  two  model  representations^  whilst  retaining  the  current 
trim  state. 

Just  as  there  are  two  methods  for  computing  the  system  response^  there  are  also  two  methods 
for  finding  a  solution  for  the  balanced  and  trimmed  state  of  the  hehcopter-load  system.  If  the 
Optimization  Toolbox  [17]  is  present  within  MATLAB®,  a  more  efficient  nonlinear  least 
squares  routine  {Isqnonlin)  is  used.  Otherwise^  the  basic  simplex  routine  included  with  the  core 
distribution  {fminsearch)  is  used.  The  nonlinear  least  squares  algorithm  iteratively  determines 
a  solution  based  on  the  accelerations^  or  external  forces  and  moments  in  each  axis.  The 
simplex  algorithm^  on  the  other  hand^  iteratively  determines  a  solution  based  on  the  sum  of 
squares  of  aU  accelerations  in  the  models  which  is  why  it  is  much  less  efficient.  For  this  reason^ 
it  is  highly  recommended  that  the  Optimization  Toolbox  be  activated. 

The  last  option  within  the  main  interface  sets  the  order  of  the  integration  algorithm  within  the 
Simulation  routine.  This  algorithm  is  a  fixed  step-size  Rimge-Kutta  integrator  and  can  be  of 
2nd/ 3rd  order  (otfe23)^  4th/5th  order  {ode45)  or  7th/8th  order  {ode78).  Obviously^  the  higher 
order  integrators  are  more  accurate;  however^  they  are  also  more  costly  in  terms  of  processing 
time  since  they  make  more  calls  to  the  force  determination  fimction.  The  lower  order 
integrators  can  also  cause  numerical  divergence  due  to  error^  particularly  for  stiff  elastic 
models^  so  a  balance  must  be  decided  upon.  The  default  order  for  new  simulation  models  is 
4th/5th^  which  has  proven  adequate  for  most  of  the  cases  tested. 

3.3  Configuration  Interfaces 

There  are  four  basic  configuration  interfaces^  which  must  all  be  set  prior  to  simulation.  These 
consist  of  the  Helicopter  Configuration,  the  Load  Configuration,  Controls  Configuration,  and  Trim 
Initialisation, 

3.3.1  Helicopter  Configuration 

The  Hehcopter  Configuration  panels  shown  in  Figure  3.3^  incorporates  options  to  select  the 
Helicopter  type,  the  Simulation  model  to  be  used^  and  the  corresponding  Weight  and  Inertia, 
Supported  Helicopter  types  include  the  BeU  UH-IH  Iroquois^  Boeing  CH-47D  Chinook^  and 
Sikorsky  CH-53D  Sea  Stallion.  Both  the  UH-IH  and  the  CH-53D  have  one  hook  and  the 
CH-47D  has  three  hooks  for  the  carriage  of  external  loads. 
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GHSLSIM:  Helicopter  Configuration 


Helicopter  type 


Simulation  model 


Weight  (lb) 


Inertia  (slugs/ff^3) 


•  UH-1H 
0  CH-47 
m  CH-53 

m  HSL 
0  HEX 


W.y. 


-Inl  xi 


I  Stability  ^  Control  Augmentation  System  (SCAS) 


lyy 


Izz 


Ikz 


F igure  3,3:  irfe / icop  ter  Configura  t ion  In te rface 


The  Simulation  model  can  be  selected  from  a  simple  Linearised  state-space  model  (HSL)^  or  a 
full  nonlinear  models  based  on  the  ROTORGEN  core  (MEX).  As  mentioned  in  Section  2.2^  the 
linearised  helicopter  model  is  only  an  approximation  and^  although  considerably  faster^  it  is 
much  less  accurate  than  the  full  nonlinear  model.  The  MEX  model  also  incorporates  aStabil  ity 
&  Control  Augmentation  System^  which  can  be  switched  on  or  off  via  a  checkbox  on  the  panel 

Text  fields  for  the  V\/eight  and  Inertia  also  allow  for  adjustment  of  these  attributes  from  within 
the  graphical  interface.  The  inertias  include  primary  axes  components  Ixxr  Iyy^  Izzr  and  Ixz-  If 
desired/  it  is  possible  to  set  any  of  these  to  infinity  (inf)/  thereby  constraining  rotation  in  the 
corresponding  axis.  This  will  be  explained  in  more  detail  in  Section  3.7.2. 

3.3.2  Load(s)  Configiu'ation 

Up  to  three  loads  can  be  configured  via  the  load  interface/  shown  in  Figure  3.4.  For  each  load/ 
there  are  options  for  the  Load  type  and  Sling  configuration^  as  well  as  fields  for  the  Load  rjjeight^ 
Load  inertia^  Size^  Attachment  locations.  There  are  also  fields  for  the  Cable  lengthy  Stiffness  and 
Damping.  Load  types  currently  supported  are  boX/  cylinder  and  plate  geometries.  The  box  load 
has  the  most  accurate  aerod5mamic  data/  based  on  a  combination  of  experimental  and 
theoretical  data/  whereas  the  cylinder  and  plate  loads  are  based  on  theoretical  data  only. 

The  Load  zueight  and  Load  inertia  in  each  primary  axis  can  be  edited  in  the  corresponding  text 
boxes  and/  as  with  the  helicopter/  they  may  be  set  to  infinity  in  order  to  constrain  the  rotation. 
AlsO/  the  size  of  the  load  can  be  specified/  although  this  is  merely  for  display  purposes  and 
has  no  effect  on  the  resulting  dynamics  of  the  system. 
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Figure  3.4:  Load(s)  Configuration  Interface 


For  the  Sling  configuration,  the  user  has  a  choice  of  four  basic  configurations,  illustrated  in 
Figure  3.5.  These  are:  Single  point.  Multiple  point,  Bifilar  and  Tandem.  Each  configuration  has  a 
different  number  of  (load)  attachments,  the  offsets  for  which  are  listed  next  to  the 
configuration  buttons.  The  Single  point  configuration  has  one  attachment,  the  Multiple  point 
has  four,  the  Bifilar  has  two  and  the  Tandem  has  four  as  well.  It  is  crucial  to  note  that  the  order 
in  which  the  attachments  are  listed  is  important.  For  two  attachment  points,  the  forward 
attachment  must  be  listed  first  and  the  aft  second.  For  four  attachment  points,  the  order  is: 
forward-starboard,  aft-starboard,  aft-port,  forward-port.  Each  offset  is  specified  relative  to  the 
load  axes,  whose  origin  is  located  at  the  load  eg,  so  most  standard  configurations  (with 
attachments  above  the  load  eg)  will  have  equal  negative  z-offset  values. 


Tandem 


Figure  3.5:  Sling  Configurations 
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With  the  graphical  interface,  the  corresponding  attachment  points  on  the  hehcopter  are 
automatically  determined.  Table  3.1  outlines  the  hehcopter  attachment  point  index  for  various 
numbers  of  loads  and  each  load^s  shng  configuration.  The  hehcopter  attachment  points,  or 
hooks,  are  numbered  from  one  to  three  in  the  order:  forward,  middle,  aft.  So,  for  example,  two 
bifilar  loads  would  occupy  hooks  (1  &  2)  and  (2  &  3).  One  tandem  load  would  occupy  hooks 
(1  &  3).  One  single  point,  one  bifilar,  and  one  multiple  point  would  occupy  hooks  1,  (2  &  3) 
and  3  respectively,  and  so  on.  As  mentioned  previously,  it  is  possible  to  define  the  attachment 
points  exphdtly  in  a  script  to  construct  the  model. 


Table  3.  T.  Helicopter  Attachment  Point  (Hook)  Index 


Load 

Configuration 

Number  of  loads  | 

1 

2 

3 

Single  point 

1 

1,3 

1,2,3 

Multiple  point 

1 

1,3 

1,2,3 

Bifilar 

1,3 

(1,2),  (2, 3) 

(1,  2),  (2, 3),  (2,  3) 

Tandem 

1,3 

(1,  2),  (2, 3) 

(1,2),  (2, 3),  (2,  3) 

As  well  as  the  sling  configuration,  the  cable  lengths  and  their  elastic  properties  must  be  set. 
Using  the  graphical  interface,  the  Cable  lengthy  Stijfness  and  Damping  for  each  cable  are  set  to 
the  same  values  specified.  Once  again,  it  is  possible  to  set  these  individually,  as  wdl  be 
explained  in  Section  3.7.2.  The  last  checkbox  shown  in  the  load  configuration  interface  is  to 
enable  Cable  slacks  which  imposes  a  zero-force  when  a  cable  is  imder  contraction.  Normally,  all 
cables  will  be  under  tension,  so  this  parameter  will  only  tend  to  become  important  when  the 
load  has  become  ims table  and  diverged.  Load  boimce  is  also  possible,  but  this  is  usually  the 
result  of  phot-induced  oschlations,  which  do  not  occur  in  the  simihation  model.  However,  the 
main  reason  that  this  option  has  been  left  imchecked  as  a  default  is  because  the  balance  and 
trim  optimisation  routine  can  have  difficulty  in  finding  a  solution  for  cables  with  slack.  This  is 
due  to  the  discontinuity  at  the  point  of  slack  and  invariable  nature  of  the  cable  force  when 
contracted.  If  every  cable  becomes  slack  (zero  cable- force)  during  balance  or  trim,  then  the 
optimisation  direction  becomes  indeterminable  and  the  routine  whl  falter. 

3.3.3  Controls  Configuration 

The  Controls  Configuration  panel  (Figure  3.6)  is  shghtly  different  than  the  previous  two 
configuration  panels,  in  that  it  ahows  the  user  to  construct  control  input  sequences  by  adding 
any  number  of  basic  input  types. 
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GHSLSIM:  Controls  Configuration 


Figure  3.  6\  Controls  Configuration  Interface 


There  are  four  tabs  -  one  for  each  control  device  -  including  Longitudinal  stick,  Lateral  stick, 
Pedals  and  Collective  Stick.  The  In-put  Types  can  be  a  step,  pulse,  doublet,  rate-feedback  or  none  and 

each  input  has  a  set  of  parameters  to  define  its  shape.  These  can  be  explained  more  clearly 
using  the  diagrams  in  Figure  3.7.  A  step  input  is  the  simplest  and  defined  by  a  Start  Time, 
Ramp  time  and  LAagnitude.  These  are,  respectively,  the  time  at  which  the  input  is  started,  the 
time  taken  to  reach  maximum  amphtude  and  the  value  of  that  maximum  amplitude.  A  pulse 
input  has  those  same  three  parameters,  plus  a  Liold  time,  which  defines  the  length  of  time  the 
control  is  held  at  its  maximum  position  (magnitude).  A  doublet  is  essentially  two  pulse  inputs 
-  one  positive  pulse  immediately  followed  by  a  negative  one  -  and  the  Start  Time,  P.amp  time 
and  iClagnitude  are  all  the  same  as  previously  defined. 


Step  ramp  time 


start  time 


magnitude 
hold  time  ramp  time 


Rate-feedback 


hold  time 

< - ► 


2  *  ramp  time  ramp  time 
hold  time 


2  *  magnitude 


Figure  3.  7.  Control  Input  Types 
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The  rate-feedback  input  type  is  quite  different  to  the  former  three  inputs.  It  is  defined  by  a  Start 
time^  a  Hold  time  and  a  Rate- gain  factor.  The  Hold  time  basically  defines  the  period  during 
which  the  rate-feedback  control  is  active  and  the  Rate  gain  defines  the  magnitude  of  control 
feedback  employed. 

This  input  is  a  simple  automatic  control  that  augments  the  commanded  controls  defined  using 
the  above  inputs.  Equation  2  describes  the  algorithm  used: 

Si,  =K^(q-0.5^da) 

Sa=K2ip-0.5*d^) 

S,=K2(r) 

S,=K,i^ 

where  Si,  denotes  the  longitudinal  stick  control^  denotes  the  lateral  sticky  S^.  denotes  the 
pedals  and  denotes  the  collective  stick.  Ki^  . .  .^  JQ  are  the  rate-feedback  gains  for  each 
control.  The  angular  body  rates  in  pitchy  roU  and  yaw  are  and  r^  respectively  and  the  rate 
of  descent  is  Si,  The  pitch  and  roU  angles  perturbed  from  their  trim  states  are  denoted  by  ^6* 
and  .  This  rate-feedback  law  has  proven  sufficient  for  steady-level  flight  application^  such 
as  initial  correction  of  trim  and  recovery  from  a  commanded  manoeuvre. 

Each  control  input  sequence  is  constructed  by  repeatedly  adding  basic  input  commands  with 
the  Acid  to  Control  button  and^  if  necessary^  clearing  the  entire  input  sequence  with  the  Clear 
Control  button.  The  basic  inputs  are  superimposed  on  to  the  current  control  sequence  and  any 
periods  of  rate-feedback  are  marked  (shaded).  Since  the  rate-feedback  essentially  replaces  any 
commanded  input  when  it  is  active^  it  is  unnecessary  to  have  any  such  input  during  those 
periods.  However^  the  rate-feedback  can  be  utilised  for  manoeuvre  recovery.  So^  for  example^ 
rather  than  employing  a  pulse  inpub  one  could  simply  use  a  step  input  followed  by  a  rate- 
feedback  inpub  Iho  automatic  control  will  tend  to  return  the  system  and  controls  to  a  state 
close  to  their  trim  state  anyway. 

Note  that  the  actual  commanded  input  or  rate- feedback  for  each  control  is  offset  by  its  trim 
position^  as  determined  during  the  trim  stage.  This  must  be  kept  in  mind  when  generating  the 
control  input  sequences^  since  the  controls  may  exceed  their  limits  when  adjusted  by  their 
trim  values  and  be  trimcated  accordingly. 

3.3.4  Trim  Initialisation 

The  last  configuration  interface  sets  fhe  initial  conditions  prior  to  trim.  Shown  in  Figure  3.8^ 
these  variables  set  the  "initial  guessb  which  the  balance  and/or  trim  optimisation  stages  begin 
from.  Fields  include:  the  Initial  Attitude,  incorporating  bank,  pitch  and  azimuth  angles;  and  the 
Initial  Controls,  Longitudinal  stick,  Lateral  stick.  Pedals  and  Collective  stick.  The  default  for  all  of 
these  is  zero,  except  the  Collective  stick,  which  is  set  to  half  of  its  upper  limit. 
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Figure  3. 8:  Initial  Conditions  Interface 


3*4  Simulation  Controls 

Oiice  the  hehcopter  and  loads  have  been  specified^  and  the  initial  conditions  set  the 
hehcopter-load  systemcanbe  Balanced  and  Trimmed.  Then^  assuming  that  the  control  inputs 
have  been  constructed^  the  simulation  may  be  Run.  Strictly  speakings  itis  not  necessary  for  the 
s  ystem  to  be  balanced^  or  trimmed  before  running.  Howev  er^  it  is  strongly  advis  ed  that  thes  e 
stages  be  followed  in  order^  to  niinirnise  disturbance  during  the  initial  period  of  simulation 
The  three  stages  under  the  heading  of  Simulation  are  described  in  the  following  section 

At  any  point  during  the  balance/ trim/ simulation  stages^  the  system  state  and  output  if 
present  can  be  saved  using  the  file  menu.  Infact^  it  is  recommended  that  the  user  saves  to  (a 
different)  file  regularly^  should  they  wish  to  repeat  a  certain  stage  without  having  to  re-run  all 
preceding  stages. 

3.4.1  Balance  Loads 

Because  the  trim  optimisation  can  be  a  difficult  task  for  multi-body  systems^  itis  often  wise  to 
balance  the  load(s)  beneath  the  helicopter  before  trirnming  the  complete  configuration 
Without  doing  so^  the  optimisation  routine  may  get  trapped  in  local  nunirna^  or  take  an 
excessively  long  time  to  reach  convergence. 

Essentially^  this  entails  holding  the  hehcopter  to  some  fixed  position^  and  attitude  while 
conducting  an  external  force  niinirnisation  on  each  load  until  they  reach  equilibrium  The  load 
subs  ystem  is  significantly  easier  to  balance^  and  when  converged^  provides  an  initial  guess  of 
the  system  state  that  is  generally  much  closer  than  previously  set  It  is  also  possible  for  the 
user  to  switch  between  inelastic  and  elastic  sling  models^  thereby  simplifying  the  shng-load 
model  even  further  for  initial  balance  optimisation 

Figure  3.9  shows  the  REPLAY  display  during  the  load  balance  stage.  The  Starts  Pause  ^  and  Stoj) 
buttons  are  all  deactivated  during  this  iterative  procedure^  and  the  display  frequency  slider 
ifreqn)  and  Tip  Trails  checkbox  have  no  effect  here.  However^  the  viewpoint  azimuth  (azim) 
and  elevation  (eleu)  shders^  and  the  Erasemode  checkbox  can  all  be  adjusted  at  anytime.  Along 
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with  the  main  three  windows  (Display,  Menu  and  Controls),  a  window  displaying  trace  plots  of 
selected  variables  that  are  used  in  the  optimisation  is  shown.  For  balance,  these  include  the 
load  angles,  ^ ,  0 ,  and  ip ,  as  well  as  the  logarithm  of  the  out-of-balance  force  component,  F. 
This  last  variable  gives  an  indication  of  the  system's  balance  during  optimisation  and  should 
converge  to  zero.  The  balance  optimisation  can  be  stopped  at  any  time  by  simply  pressing  the 
STOP  button  in  the  variable-tracking  window. 


Figure  3.9:  Loads  Balance  Display 


If  the  Verbose  Output  flag  has  been  checked  in  the  program  preferences,  the  optimisation 
routine  writes  information  to  the  MATLAB®  command  window  as  shown  in  Figure  3.10.  For 
the  balance  and  trim  optimisation,  this  information  includes  the  iteration  number,  the 
cumulative  number  of  function  evaluations,  the  residual,  step  size  and  directional  derivative. 
The  residual  should  approach  zero,  and  both  the  step  size  and  directional  derivative  should 
get  smaller  as  the  iterations  progress.  More  detail  on  these  indicators  is  available  in  the 
Optimization  Toolbox  documentation  [17].  When  the  optimisation  has  finished,  the  final 
values  for  each  of  the  variables  shown  in  the  variable  tracking  window  are  listed. 
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Figure  3.10\  MATLAB®  Command  Window  Output 

A  suggested  procedure  for  balancing  the  load(s)  under  the  helicopter  is  as  follows: 

1.  First  set  the  initial  conditions  for  the  helicopter,  if  a  good  estimate  is  known,  using  the 
Trim  Initialisation  interface. 

2.  Set  the  Sling  model  to  inelastic  in  the  Main  interface,  then  balance. 

3.  If  the  load  has  balanced  adequately,  set  the  Sling  model  to  elastic  and  balance  again. 
Otherwise,  re-run  the  balance  optimisation,  possibly  modifying  the  initial  conditions 
beforehand.  This  may  be  necessaiy  if  the  routine  is  having  trouble  converging  to  a 
balanced  state. 

If  one  does  not  require  an  elastic  simulation,  then  the  third  step  is  imnecessaiy.  Once 
balanced,  the  user  may  go  on  to  trim  the  complete  system.  Since  the  balance  stage  operates 
solely  on  the  load  and  does  not  call  the  ROTORGEN  helicopter  model,  it  is  less  prone  to 
accuracy  limitations.  However,  these  may  be  encoimtered  dining  the  trim  stage,  as  will  be 
explained  next. 

3.42  Trim  System 

Although  not  necessaiy,  it  is  advisable  to  trim  the  helicopter-load  system  prior  to  simulation 
to  avoid  any  initial  disturbance  due  to  out-of-balance  forces.  Depending  on  the  fidelity  of  the 
simulation  required,  it  is  possible  to  utilise  an  inelastic  sling  model,  or  an  elastic  one,  and  the 
trim  routine  can  handle  either.  The  REPLAY  display,  shown  in  Figure  3.11,  is  much  the  same 
as  that  rendered  in  the  balance  stage,  apart  from  the  addition  of  the  control  input  variables. 
These  include,  in  order,  the  longitudinal  stick  <5^ ,  the  lateral  stick  ,  the  pedals  ,  and  the 
collective  stick  .  The  controls  window  is  also  activated  during  this  stage,  for  a  clearer 
indication  of  the  control  stick  displacements.  They  are  generic  slider-type  controls  and 
incorporate:  the  longitudinal  and  lateral  stick,  placed  vertically  and  horizontally  in  the  main 
axes;  the  pedal  controls,  placed  horizontally  below  the  main  axes;  and  the  collective  stick, 
placed  vertically  to  the  right  of  the  main  axes.  The  controls  are  normalised  with  respect  to 
their  hard  limits  in  aU  of  the  display  axes. 
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Figure  3.11:  System  Trim  Display 


If  desired,  one  can  change  the  viewpoint  azimuth  and  elevation,  and  the  Erasemode  during 
trim  using  the  controls  in  the  Menu,  and  as  with  the  balance  stage,  it  is  possible  to  stop  the 
trim  optimisation  at  any  time  using  the  button  displayed. 

The  suggested  procedure  for  trimming  the  helicopter-load  system  is: 

1.  First  set  the  initial  conditions  for  the  helicopter,  and  balance  the  load  under  it. 

2.  Trim  the  system  using  either  in  inelastic  sling  model,  or  an  elastic  one  as  required.  It  is 
not  necessaiy  to  do  both. 

3.  If  the  system  has  not  reached  trim,  re-run  the  trim  optimisation.  Check  that  the  load  is 
balanced.  If  the  system  still  does  not  trim,  try  a  different  set  of  initial  conditions. 
Failing  this,  it  may  be  necessaiy  to  fly  the  helicopter  and  load  into  an  approximate 
trim  state  through  simulation.  See  Section  3.7.2  for  a  further  discussion. 
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Note  that  for  some  configurations^  it  may  be  difficult  to  find  a  satisfactory  trim  point.  For 
example^  a  lightweight  load  flown  imder  a  hehcopter  at  high  speed  will  tend  to  be  very 
ims table.  Moreover^  it  may  not  be  physically  possible  to  attain  trim  at  all^  such  as  the  case 
when  a  load  is  too  heavy  for  the  hehcopter  to  hft.  The  user  should  be  aware  of  the  limitations 
specific  to  each  configuration  before  attempting  trim  or  simulation  using  GHSLSIM. 

Unfortunately^  there  are  also  accuracy  limitations  in  the  ROTDRGEN  (MEX)  hehcopter  model. 
MATLAB'®^  and  the  models  generated  within  the  environment  such  as  the  loads^  are  all 
intrinsically  double-precision^  whereas  the  ROTORGEN  model  is  much  less  accurate.  As  a 
consequence^  truncation  error  may  cause  the  trim  optimisation  to  get  ^stuck^  in  a  discrete  state^ 
or  it  may  not  converge  to  a  satisfactory  solution.  Often  however^  the  trim  optimisation  whl 
reach  a  satisfactory  state  within  the  maximum  number  of  iterations^  even  though  the  function 
may  suggest  the  optimisation  be  re-rtm  (possibly  increasing  the  maximum  number  of  function 
evaluations).  If  it  appears  to  have  done  so,  according  to  the  residual  error,  then  it  can  be 
considered  correctly  trimmed. 

3.4.3  Run  Simulation 

The  REPLAY  window  appears  again  for  the  simulation  stage,  shown  in  Figure  3.12.  In  the 
variable  tracking  window,  the  hehcopter  inertial  velocities  ( j^S,  3S,  &)  are  plotted,  as  weh  as  the 
hehcopter  angles  and  control  inputs.  The  REPLAY  function  is  still  cahed  iteratively,  so  none  of 
the  main  menu  commands  function,  however,  a  counter  window  provides  an  indication  of  fhe 
simulation  progress  and  the  simulation  can  be  stopped  at  any  time  as  with  the  trim  stage.  In 
addition,  if  the  Verbose  Output  flag  has  been  checked  in  the  program  preferences,  the 
integration  routine  will  output  the  time,  the  absolute  hehcopter  velocity,  the  relative  load 
velocity  and  an  estimate  of  the  integration  error  in  the  MATLAB®  command  window. 

Typicahy,  the  system  is  too  complex  to  rtm  in  real-time  and  therefore,  simulations  with  a 
small  times tep,  or  a  high  order  integration  algorithm  wtil  take  longer  to  complete.  However,  if 
the  simulation  has  diverged  due  to  trimcation  error,  it  may  be  necessary  to  increase  the  order, 
or  reduce  the  times  tep. 
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Figure  3.12:  Simulation  Run  Display 

3.5  Output  Controls 

The  output  controls  are  intended  to  be  used  following  simulation,  and  include  buttons  to 
Replay  Simulation,  Plot  Time-Histories,  and  Generate  VRML  code.  Both  the  simulation  replay  and 
VRML  model  provide  a  real-time  three-dimensional  graphical  display  of  the  simulation, 
whereas  the  time-history  plots  are  intended  for  analysis  of  the  output  data. 

3.5.1  Simulation  Replay 

This  fadhty  opens  the  same  REPLAY  windows  as  for  the  simulation  stages;  however,  this 
time  it  is  called  with  the  full  position,  orientation  and  control  data  sets,  rather  than  just  the 
data  at  any  particular  step  in  an  iterative  fashion.  Executing  the  REPLAY  function  in  this  way 
allows  full  control  of  the  Start,  Pause,  and  Stop  cornmands  in  the  menu,  as  well  as  several 
options  not  previously  available.  An  example  of  the  display  is  shown  in  Figure  3.13. 
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Figure  3.13:  Simulatkm  Replm/ 


Freqn  denotes  the  integer  display  frequency  for  the  simulation,  that  is,  the  number  of  frames 
skipped  before  each  display  step.  For  example,  a.freifn  of  1  will  display  every  frame,  and  a 
frecpt  of  10  will  display  every  10th.  The  viewpoint  can  also  be  moved  about  360^  in  each  axis 
using  the  azim  andelev  sliders.  If  the  Erasemode  option  is  checked  (the  default),  each  frame  will 
be  erased  before  rendering  the  next.  If  it  is  unchecked,  a  trail  or  history  of  frames,  such  as  that 
shown  in  the  figure  will  be  displayed.  Note  that  these  frames  aren't  permanent  objects,  and 
once  the  viewpoint  changes  the  trail  will  disappear.  Another  option  to  help  with  the 
visualisation  is  to  display  Tip  Trails  during  the  simulation  replay.  The  Tip  Trails  are  simply 
coloured  lines  that  follow  the  port  and  starboard  extremities  of  the  helicopter,  much  like  long 
streamers.  Below  the  controls  and  options  is  a  frame  displaying  the  current  time  and  timestep. 
The  associated  slider  can  be  moved  to  any  point  within  the  simulation  range  and  the 
helicopter-load  model  will  be  set  to  the  corresponding  position.  The  last  box  in  the  menu 
window  is  the  Realtime  indicator.  During  replay,  regardless  of  the  display  frequency,  this 
changes  colour,  depending  on  whether  the  simulation  is  being  replayed  in  real  time:  green 
indicates  real  time  and  red  indicates  a  refresh  rate  slower  than  real  time. 

3.5.2  Time-History  Plots 

The  time-history  plots  give  a  qualitative  overview  of  the  simulation  run  for  the  helicopter  and 
load(s).  The  variables  that  are  plotted  depends  on  the  Simulation  axes  defined  in  the  simulation 
configuration,  as  well  as  the  control  inputs  and  the  number  of  loads.  An  example  of  one  set  of 
plots  is  shown  in  Figure  3.14,  for  a  single- load  configuration. 
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Figure  3. 14\  Time-History  Plots 

Here,  the  body-axes  velocities  (w,  v,  w)  and  the  body-axes  angular  rates  (p,  q,  r)  are  included, 
as  well  as  the  Euler  angles  ( ^ ^  ^  )  and  the  four  control  inputs  5^,  5^,  5^  ). 
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3.5.3  VRML  Generation 

The  Virtual  Reality  Modelling  Language  (VRML)  is  an  ideal  language  for  the  visualisation  of 
dynamic  flight  simulations  [18].  Using  the  VRML  Generation  command,  the  user  can  export 
the  simulation  in  VRML  format,  which  can  then  be  viewed  and  replayed  in  a  VRML  browser. 
Figure  3.15  is  a  screen  capture  of  the  Internet  Explorer  web  browser  with  the  Cortona  VRML 
plug-in  installed. 


Figure  3.15:  VRML  Simulation  Replay 


It  is  important  to  make  sure  the  VRA4L  Directory  is  set  correctly  in  the  program  preferences, 
that  is,  the  directory  in  which  all  the  VRML  helicopter  and  load  models  reside.  Otherwise, 
they  will  not  appear  in  the  VRML  simulation. 

To  run  the  simulation  replay,  the  user  simply  clicks  on  either  the  helicopter  or  load.  It  is  not 
possible  to  pause  or  stop  the  simulation  midway  through  the  run,  and  to  go  back  to  the  first 
frame  without  restarting,  the  user  must  reload  the  model  (using  the  browser's  reload 
command). 

There  are  several  viewpoints  to  choose  from,  including  User,  Side,  Front,  Top,  and  Behind. 
These  are  all  static,  in  the  sense  that  their  inertial  position,  relative  to  the  helicopter,  is  fixed 
throughout  the  simulation,  whilst  the  view  angle  remains  constant.  The  dynamic  viewpoints 
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include  Fixed-Inertial,  Fixed-Zoom-Pan,  Fixed-Aircraft  and  Fixed-Down,  Unlike  the  first  set,  the 
view  angle  for  these  changes  during  the  simulation.  The  Fixed-Inertial  viewpoint  represents 
one  whose  position  is  fixed  in  space,  such  as  an  observer  on  the  ground.  Fixed-Zoom-Pan  is  the 
same,  except  the  focal  length  changes  with  time  such  that  the  helicopter-load  configuration 
occupies  the  full  view  throughout.  The  Fixed-Aircraft  and  Fixed-Down  viewpoints  have  fixed 
offsets  from  the  helicopter,  which  remain  constant  and  their  view  angles  change  with  any 
variation  in  the  helicopter  angle.  The  former  is  located  at  a  position  and  angle  similar  to  the 
User  viewpoint,  while  the  latter  is  located  just  below  the  middle  helicopter  attachment  point, 
looking  down  at  the  load.  As  well  as  the  listed  viewpoints,  it  is  possible  to  navigate  the 
viewpoint  to  any  position  and  angle  in  space  using  the  VRML  browseT  s  controls.  When  doing 
so,  it  is  important  to  remember  that  the  variation  of  the  viewpoint  position  and /  or  angle  is  not 
altered,  just  the  offset.  So,  for  example,  changing  the  initial  position  of  the  User  viewpoint  will 
follow  the  helicopter  as  expected,  but  changing  the  initial  position  of  the  Fixed-Inertial 
viewpoint  may  result  in  a  strange  replay  that  does  not  follow  the  helicopter  at  all. 

In  addition  to  the  helicopter,  load  and  cables,  a  load  displacement  patch  is  superimposed  on 
the  VRML  model  to  aid  with  visualisation  of  the  load's  swinging  motion.  This  patch  extends 
from  the  helicopter  attachment  point  to  the  load  centre-of- gravity  and  indicates  the  amoimt  of 
relative  displacement  seen  by  the  load  from  its  trim  position. 

3.6  Example:  3000  L  Water  Tank 

An  example  is  given  here  for  the  CH-47D  and  a  3000  L  water  tank,  suspended  by  four  slings, 
as  shown  in  Figure  3.16.  This  is  a  simple  load  and  can  be  adequately  represented  by  a  box¬ 
shaped  container.  The  centre  of  gravity  is  assumed  to  remain  at  the  centre  of  the  tank  at  alt 
times. 


. . . miimiiiLi*-.. 
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POTABLE  WATER  ONLY 


Figure  3.16:  3000  L  Water  Tank 
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The  control  inputs  specified  include  a  rate-feedback  on  the  longitudinal  stick  to  keep  the  nose 
level  throughout  the  manoeuvre,  as  well  as  a  doublet  on  the  lateral  stick  and  a  pulse  on  the 
pedals  to  excite  the  load.  For  this  case,  they  are  not  intended  to  accurately  represent  any  real- 
world  inputs  that  might  typically  be  employed,  but  specifically  designed  to  excite  the  load 
imder  the  heUcopter. 

Below  is  a  list  of  the  steps  to  configure,  trim  and  simulate  the  hebcopter  load  system.  They 
represent  just  one  approach,  especially  for  trim,  that  can  be  used  to  achieve  the  same  outcome. 
For  example,  another  scheme  to  obtain  trim  would  be  to  first  balance  the  inelastic 
configuration,  then  trim  the  inelastic  configuration,  and  then  trim  the  elastic  configuration. 

The  steps  are  as  follows: 

1.  Start  the  program,  by  typing  at  the  command  prompt: 

»  ghs1s1m 

2.  First,  check  the  program  preferences.  Select  File  -►  Preferences. . .  from  the  main  menu. 
Confirm  that  the  Web  browser,  VRML  Directory  and  GZIP  Executable  fields  are  correct 
and  if  not,  update  them  accordingly.  Check  all  of  the  checkboxes  for  a  detailed  output 
and  visualisation. 

3.  Next,  in  the  main  interface,  set  the  following  fields: 


Simulation  time  (s) 

10 

Integration  timestep  (s) 

0.1 

Number  of  steps 

101 

Trim  airspeed  (ft/s):  Forward 

100 

Trim  altitude  (ft) 

100 

Simulation  axes 

(s)  combined 

Angular  representation 

®  Euler 

Simulation  type 

®  nonlinear 

Sling  model 

(s)  inelastic 

Trim  routine 

®  Isqnonlin 

Simulation  routine 

®  ode45 

This  will  create  a  10  second  simulation  with  101  timesteps.  The  forward  speed  is 
approximately  60  kn  and  the  altitude  100  ft.  AU  options  are  set  to  their  defaults,  except 
for  the  Sling  model,  which  is  set  to  inelastic  for  the  balance  stage. 

4.  Press  the  Configure  helicopter  button.  In  the  Hehcopter  Configuration  window  that 
pops  up,  type  in  the  f oho  wing  options: 

Helicopter  type  ®  CH-47 

Simulation  model  ®  MEX 

For  this  rim,  leave  the  Stability  &  Control  Augmentation  System  (SCAS)  checkbox 
imchecked. 

Weight  (lb)  46000 

Inertia  (slug.ft^2):  Ixx;  lyy;  Izz  43200;  250000;  234000 

Press  the  OK  button  to  return  to  the  main  interface.  These  are  the  hgures  for  a  fuhy 
loaded  Chinook  CH-47D,  with  the  SCAS  switched  off. 
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5.  Press  the  Configure  loads  button.  In  the  Hehcopter  Configuration  window  that  pops  up, 
type  in  the  following  options  and  fields: 


Load  type 

Load  weight  (lb) 

Load  inertia  (slug.ft^2):  Ixx;  lyy;  Izz 

Size  (ft):  Length;  Width;  Height 

Sling  configuration 

Cable  length  (ft) 

Stiffness  (Ibf/ft) 

Damping  (Ibf/ft/s) 

Attachments  (ft):  x-offset;  y-offset;  z-offset 

®  box 

3330 

1000; 4000;  4000 

8.00;  6.46;  3.84 
(s)  Multiple  point 

19.5 

20000 

500 

4.00;  3.23; -1.92 
-4.00;  3.23; -1.92 
-4.00; -3.23; -1.92 

4.00;  ^.23; -1.92 

Leave  the  Cables  slack  (zero  force)  under  contraction  checkbox  unchecked.  This  represents 
an  empty  3000  L  water  tank,  supported  by  four  shngs  attached  to  each  upper  corner  of 

the  load.  Since  there  is  only  one  load  with  a  multiple 

-point  configuration,  the  shngs 

will  automatically  be  located  to  the  hehcopter  middle  cargo  hook.  At  this  stage,  there 
is  no  attempt  to  implement  the  correct  centre-of-gravity  position,  and  consequently,  it 
is  placed  at  the  centre  of  the  load.  Press  the  OK  button  to  return  to  the  main  interface 

again. 

Press  the  Configure  controls  button.  In  the  Controls  Configuration  window  that  appears, 
make  sure  the  Longitudinal  stick  tab  is  active  and  type  in  the  following  options: 

Input  type 

(s)  rate-feedback 

Start  time  (s) 

0 

Hold  time  (s) 

10 

Rate  gain 

40 

Now  press  the  Add  to  Control  button  to  add  the  specified  input  to  the  longitudinal 
control.  The  entire  control  graph  should  fill  with  grey,  indicating  the  region  of  control 
feedback.  Next,  press  the  Lateral  stick  tab  to  activate  that  control  and  type  in  the 

foUowing: 

Input  type 

(§)  doublet 

Start  time  (s) 

1 

Ramp  time  (s) 

0.5 

Hold  time  (s) 

1 

Magnitude  (in) 

0.2 

Press  the  Add  to  Control  button  to  add  the  specified  input  to  the  lateral  control.  A 
doublet  input  with  those  properties  wiU  be  drawn  to  the  control  graph.  Then  within 

the  same  frame  (Lateral  stick\  type  in  the  following: 

Input  type 

®  pulse 

Start  time  (s) 

4.5 

Ramp  time  (s) 

0.5 

Hold  time  (s) 

0.5 

Magnitude  (in) 

-0.2 
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Again,  press  the  Add  to  Control  button  to  superimpose  the  specified  input  to  the  lateral 
control.  Adding  this  pulse  effectively  lengthens  the  second,  negative  component  of  the 
doublet-type  control  sequence.  Next,  press  the  Pedals  tab  to  activate  that  control  and 
type  in  the  following: 


Input  type  (§)  pulse 

Start  time  (s)  1 

Ramp  time  (s)  0.5 

Hold  time  (s)  2 

Magnitude  (in)  0.1 


Press  the  Add  to  Control  button  to  add  the  specified  input  to  the  pedal  controls.  A  pulse 
input  wlU  be  drawn  to  the  graph.  Next,  press  the  Collective  stick  tab  to  activate  that 
control  and  type  in  the  following: 

Input  type  ®  none 

Press  the  Add  to  Control  button  to  add  the  null  input  to  the  collective  controls.  Lastly, 
press  the  OK  button  to  return  to  the  main  interface. 

7.  Press  the  Initial  Conditions  button.  In  the  Trim  Initiahsation  window  that  appears,  type 
in  the  following  options: 

Initial  Attitude  (deg):  bank;  pitch;  azimuth  0;  5;  0 

Initial  Controls  (in):  Longitudinal;  Lateral;  0;  0 

Pedals;  Collective  0;  4.56 

The  hehcopter  should  have  a  shght,  positive  pitch  angle  at  trim,  hence  the  initial  guess 
of  5°.  Each  of  the  controls  is  typically  set  to  its  central  (range-average)  positions.  For 
the  Collective  stick,  which  ranges  from  ff'  to  9.12"  for  the  CH-47D,  this  is  4.56",  and 
for  the  rest  it  is  zero.  When  finished,  press  the  OK  button  to  return  to  the  main 
interface. 

8.  Save  the  simulation  configuration  data  via  the  main  menu:  File  -►  Save  As.  It  is  highly 
recommended  that  the  user  saves  at  this  point,  in  order  to  avoid  having  to  re-enter  all 
of  the  above  data  should  anything  go  wrong!  This  whl  save  all  the  configuration  data 
and  the  initial  trim  state  guess.  Now,  if  the  configuration  needs  to  be  altered,  or  any 
subsequent  balance,  trim  or  simulation  re-rim,  it  is  simply  a  matter  of  loading  the 
saved  file  via  the  same  menu  (File  Open),  and  resuming  from  there. 

9.  First  step  of  the  simulation  stage:  Balance  the  loads  by  pressing  the  BALANCE  LOADS 
button.  This  first  simple  inelastic  configuration  should  balance  within  a  few  iterations. 
The  load  should  have  a  pitch  angle  of  approximately  -5°,  that  is,  swung  aft  by  5°  due 
to  the  bluff-body  drag  in  effect  at  100  ft/s. 

10.  Now  that  a  better  idea  of  the  initial  state  has  been  determined,  one  can  proceed  to  the 
elastic  configuration  with  confidence.  Set  the  Sling  model  to  elastic  with  the  radio 
button  on  the  main  interface  and  balance  again  (BALANCE  LOADS).  The  cables  wfil 
stretch  sfightly  to  coimteract  the  force  exerted  by  the  load  weight  and  drag. 
Consequently,  the  load  position  and  angle  will  change,  but  not  by  much.  The  pitch 
angle,  for  example,  will  still  be  approximately  -5°. 
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11.  Proceed  directly  to  trim,  without  making  any  further  alterations  to  the  configuration. 
Press  the  TRI M  button.  This  time,  the  hehcopter's  attitude  and  controls  will  change,  as 
well  as  the  load's  position  and  angle.  Since  the  algorithm  is  now  optimising  the  full  12 
degrees- of- freedom,  it  will  take  shghtly  longer  to  converge.  For  this  configuration,  the 
nonlinear  least  squares  trim  routine  will  take  around  50  iterations  before  the  tolerance 
criterion  is  met.  The  resulting  trim  values  displayed  for  the  hehcopter  are: 


Table  3.2:  Helicopter  Trim  State  with  3000  L  Water  Tank  Load 


phi 

theta 

psi 

del_b 

del_a 

del_r 

del_c 

-0.10 

2.95 

0.00 

1.00 

0.06 

0.16 

5.37 

Save  the  configuration  data  via  the  main  menu  as  before.  If  the  same  file  is  chosen  (or 
saved  directly,  using  File  Save),  this  will  update  the  Trim  Initiahsation  variables  to 
their  new  values.  Saving  at  this  point  can  be  considered  optional,  although  it  is 
recommended  if  multiple  simulation  runs  are  to  be  executed  from  the  same  trim  state. 
For  example,  if  the  user  is  interested  in  testing  the  effect  of  different  control  inputs  on 
the  same  configuration. 

12.  Rim  the  simulation  with  the  RUN  button.  For  the  101  time  steps,  this  takes  arotmd  30 
seconds  on  an  Intel  Pentium-4 1.9  GHz  processor  machine.  If  the  Replay  Window  option 
was  checked,  the  user  will  see  the  hehcopter  and  load  being  simulated  as  the 
integration  routine  progresses.  When  the  simulation  has  completed,  a  new  Replay 
window  will  appear  with  the  full  simulation  ready  for  replay. 

13.  Save  the  simulation  configuration  and  output  data  via  the  main  menu:  File  -►  Save  As. 
Once  again,  it  is  highly  recommended  that  the  user  saves  at  this  point. 

14.  If  the  Replay  Window  option  was  not  checked,  or  the  Replay  window  was  closed,  it  is 
still  possible  to  replay  the  simulation  using  the  Replay  Simulation  button  imder  the 
Output  controls.  Close  the  windows  using  the  button  on  the  Replay  menu  when 
finished. 

15.  Plot  the  time-histories  using  the  Plot  Time-Histories  button. 

16.  Generate  VRML  model  and  replay  with  the  Generate  VRML  button.  The  resulting 
VRML  file  will  be  automatically  saved  to  the  VRML  Directory  specified  in  the 
preferences.  Exit  the  program. 

The  configuration  can  now  be  re-trimmed,  or  the  simulation  re-rim  from  the  files  saved 
earher.  It  is  also  possible  to  replay  the  simulation  again,  or  plot  the  time-histories  from  any 
complete  rim  that  was  saved. 
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3.7  Advanced  Topics 

The  Graphical  User  Interface  (GUI)  is  sufficient  for  most  purposes^  and  provides  a  convenient 
tool  for  slung  load  simulation  and  analysis  that  requires  very  little  knowledge  of  the 
tmderlying  system.  Indeed^  the  GUI  is  just  a  higher-level  simplified  interface  to  the  full 
HSLSIM  package.  Many  options  are  hidden  and  various  assumptions  are  made  regarding  the 
configuration  details  within  the  program.  However,  it  is  possible  to  access  the  lower-level 
scripts  and  functions  of  the  package,  and  once  the  user  has  gained  some  experience  in  using  it, 
they  may  wish  to  write  fheir  own  scripts,  or  tailor  the  code  to  meet  their  requirements. 

3.7.1  Automation  &  Multiple  Simulation  Runs 

The  following  section  details  an  example  for  the  automated  execution  of  a  suite  of  simulation 
runs.  In  this  example,  the  configuration  from  the  previous  section  is  used,  that  is,  a  CH-47D 
helicopter  and  3000  L  water  tank  suspended  by  four  shngs.  The  aim  of  the  experiment  is  to 
assess  the  effect  of  airspeed  and  load  weight  on  the  stability  of  the  load  during  lateral  doublet 
manoeuvres.  To  this  end,  the  maximum  load  deviation  will  be  recorded  for  a  suite  of 
simulation  runs  over  a  range  of  trim  airspeeds  and  load  masses.  Refer  to  Appendix  A.2.1  for  a 
full  listing  of  the  code,  ch47b3k1  ts1m.  m. 

Variables  that  can  be  safely  modified  by  the  user  are  all  defined  at  the  top  of  the  file  in  the 
configuration  section.  These  include  Entitle,  BLindir,  Bl-V,  BlJI,  BlT,  Prefs,  Config  and 
wrl  f  i  1  es.  The  variables  BL-ti  tl  e  and  B_jnd1  r  specify  the  filename  prefix  and  directory  used 
when  the  data  from  each  run  is  saved  to  disk.  It  is  not  compulsory  to  incorporate  the  tide  into 
the  directory  name,  but  has  been  done  for  clearer  organisation.  The  three  vectors,  BLJH/  BlV, 
and  BlT  detail  the  flight-points  in  terms  of  load  mass,  airspeed,  and  control-input  magnitude, 
respectively,  at  which  the  configuration  is  to  be  trimmed  and  simulated.  BLJII/  consists  of  10 
points  spaced  uniformly  between  3330  lb  and  9943  lb,  BlV  consists  of  6  points  between  60  kt 
and  120  kt,  and  B^-T  consists  of  just  three  points  at  0.0, 1.0  and  2.0.  Consequently,  several  of  the 
resulting  special  variables  employed  to  store  results,  such  as  the  maximum  load  deviation 
LdA,  wdl  have  three  dimensions  and  be  of  size  BlJHN  X  BJVN  X  Bl.TN,  corresponding  to  the 
length  of  each  vector  (10  x  6  x  3).  The  Pref  S  struct  contains  general  program  preferences,  as 
specified  from  within  the  GUI.  These  include  the  VRML  Directory,  the  location  of  the  GZIP 
Executable  if  present,  and  several  display  flags.  The  Web  Browser  field  is  not  used  here,  since 
the  browser  is  not  launched  automatically  at  any  stage.  The  Conf  1  g  struct  contains  data  for 
the  simulation,  helicopter,  load,  controls  and  trim  initialisation.  Most  of  these  fields  should  be 
self-explanatory,  with  a  few  exceptions.  Firstly,  note  that  the  horizontal  trim  velocity  is 
initially  equal  to  BlV(1)  in  ft/s,  hence  the  conversion  factor.  The  load  weight  is  equal  to 
BujnCl),  and  the  load  inertia  (in  x,  y,  and  z)  is  approximated  by  a  simple  linear  function  of 
BLjn(l).  likewise,  the  magnitude  of  the  lateral  stick  and  pedal  inputs  (2nd  and  3rd  elements 
of  Ctrl  data)  are  functions  of  Bu-TC1)  .  These  are  aU updated  during  each  iteration  of  the  main 
loop.  Note  also  that  the  number  and  order  of  the  loaddata  attachment  offsets  are  important 
and  specific  to  the  type  of  sling  configuration  used.  For  a  single  attachment,  there  is  only  one 
set  of  (x,  y,  z)  locations.  For  a  bifilar  configuration,  there  are  two,  ordered  from  the  forward 
attachment  to  the  aft.  For  both  multiple  and  tandem  configurations,  there  are  four,  ordered: 
forward-starboard,  aft-  starboard,  aft-port,  and  then  forward-port  (i.e.,  clockwise  ordering). 
Both  the  helicopter  and  load  inertias  are  represented  by  the  traditional  3x3  matrices  [19]. 


28 


DSTOTN-0855 


Lastly^  pay  careful  attention  to  the  braces  used  for  each  field  -  some  are  matrices  or  scalars 
(with  square  brackets)^  and  some  are  cells  (with  curly  brackets).  The  2nd  and  3rd  elements  of 
wr  1  f  i  1  es  define  the  location  of  the  heficopter  and  load  models  to  be  referenced  in  the  output 
VRML  file.  The  1st  element  is  left  empty  in  this  section^  but  is  filled  later  with  an  automatically 
generated  filename. 

Other  variables  that  are  created  include  the  main  index  matrix  which  specifies  the  order 

of  the  trim-simulation  rims.  Since  the  initial  trim  state  guess  for  each  iteration  is  based  on  the 
previous  trim  state  determined  (through  optimisation)^  this  order  is  arranged  to  rrdnirnise  the 
difference  between  consecutive  trim  points.  Hence^  the  main  loop  first  ascends  through  the 
load  mass  sequence  at  the  rninimtim  velocity^  then  descends  through  the  load  mass  sequence 
at  the  next  highest  velocity^  then  ascends  through  the  load  mass  sequence  again  for  the  next 
highest  velocity  and  so  on.  This  zig-zag  path  over  the  mass- velocity  space  is  just  one  scheme 
that  can  be  employed  to  cover  the  range  of  flight-points.  The  main  aim  of  any  such  scheme  is 
to  start  from  the  best  possible  guess  of  trim  at  each  flight-point  and  therefore  minimise  the 
occurrence  of  imconverged  cases  and/or  errors.  The  matrix  BLdA  holds  the  maximum  load 

deviation  for  each  flight-point.  B _ holds  a  flag  for  each  flight-poinfi  where  a  value  of  0 

indicates  that  the  configuration  is  imtrimmed^  1  indicates  that  it  is  trimmed  and  2  indicates 
that  it  has  been  (trimmed  and)  simulated.  The  matrix  B_uO  holds  the  trim  configuration 
velocities^  B_rO  holds  the  trim  positions  and  orientations^  BLJCt  holds  the  trim  states  and  BLdO 
holds  the  trim  controls  for  each  flight-point. 

The  main  functions^  hsl  trim  and  hsl  si m  actually  use  a  set  of  different  variables^  as  opposed 
to  Conf  i  g  and  Pref  in  which  the  options  and  configuration  details  are  kept.  These  are  the 
global  variables  HDAT_,  LDAT^  CDAT^  opt_,  and  pref_  .  The  first  three  hold  the  equivalent 
hehcopter,  load  and  cable  data  and  the  next  two  hold  the  program  options  and  preferences. 
For  more  detail^  see  the  script  ghs_i  ni  which  converts  the  data  in  Conf  i  g  and  Pref  s  to  the 
global  variables  used. 

Upon  execution^  the  matrix  B _ is  checked  to  see  if  the  script  has  been  rim  before  and  if  nofi  it 

creates  all  of  the  above  variables  and  proceeds  to  the  main  loop.  If  it  has^  this  first  step  is 
skipped^  as  well  as  any  flight  points  encoimtered  during  each  cycle  of  the  loop  that  have  been 
successfully  trimmed  and  simulated.  If  any  particular  flight-point  has  been  trioimed^  but  not 
simulated,  the  remaining  simulation  cases  will  be  run  from  the  corresponding  trim  state.  In 
other  words,  the  program  attempts  to  trim  at  each  flight  point,  and  then,  if  successful,  rims 
each  of  the  simulation  cases  for  that  point.  For  each  flight-point  following  the  first,  the  trim 
state  used  to  initiafise  the  trim  optimisation  is  that  determined  from  the  previous  flight-point. 
If  the  system  cannot  be  trimmed  to  within  tolerance  at  any  particular  flight-point,  the  trim 
state  used  is  that  from  two  iterations  ago.  If  the  system  still  cannot  be  trimmed,  the  trim  state 
from  three  iterations  ago  is  used,  and  so  on  until  the  fist  of  all  previous  states  has  been 
exhausted.  This  effectively  reduces  the  likelihood  of  the  optimisation  being  imable  to 
converge  because  of  a  poor  initial  trim  state  guess. 
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At  every  step  of  the  main  loop^  the  load  mass  and  inertia^  and  the  trim  velocity  are  changed  to 
their  appropriate  values  as  given  by  BlhI  and  BlV.  As  mentioned  above^  the  inertia  is 
approximated  by  a  simple  linear  fimction  of  BlJH.  The  vertical  location  of  the  water  tankas 
centre- of-gravity  is  also  approximated  using  a  simple  fimction  based  on  the  water  level  for 
any  given  weight  mj .  Since  the  water  level  is  2.0  in  below  centre  for  both  empty  (mi  =  3330  lb) 
and  full  (mN  =  9943  lb),  and  the  height  of  the  tank  is  46.0  in,  the  eg  offset  can  be  expressed: 


=  25.0  - 


nii 


23.0^mi  + 


23.0  ^(nij  -  mi) 

(niN  -  mi) 


(3) 


After  setting  the  load  mass  and  inertia,  and  the  trim  velocity  for  the  current  flight-point,  the 
balance/ trim  stage  is  commenced.  This  is  done  in  three  steps:  First,  balance  the  inelastic  load 
configuration;  then  balance  the  elastic  load  configuration;  and  lastly,  trim  the  entire 
hehcopter-load  system.  For  each  balance  step,  the  following  commands  are  executed: 

1.  Shng  model  is  set  to  elastic  or  inelastic 

2.  Global  variables  are  initiahsed  with  ghs_i  nit 

3.  Load  vertical  eg  position  is  calculated  and  used  in  the  geometric  sub-struct  LDAT_.  S 

4.  Initial  load  position  is  set  directly  beneath  the  fixed  helicopter 

5.  Replay  window  is  created  if  flagged  in  the  program  preferences 

6.  Load  is  balanced  beneath  the  fixed  hehcopter  via  hsll oad 

7.  Load  positional  trim  data  is  updated 

The  trim  step  is  much  simpler,  entailing  the  following  commands: 

1.  Global  variables  are  initiahsed  with  ghs_i  ni  t 

2.  Replay  window  is  created  if  flagged 

3.  Trim  time-control  matrix  TDO  is  set  to  its  initial  value 

4.  Hehcopter  load  system  is  trimmed  viahsltrim 

5.  Trim  control  vector  dO  is  updated 

The  sum  of  squares  of  the  residual  external  force  is  then  checked  against  a  tolerance  (0.1  in 
this  example),  and  if  unsuccessful,  the  trim  index  is  decremented,  the  corresponding  element 

of  B _ is  set  to  0  and  program  control  returns  to  the  top  of  the  loop  for  another  attempt  at  trim 

starting  at  the  previous  flight-point  trim  state.  If  the  check  is  successful,  the  matrices  BLuO, 

B_rO,  BLdO  and  B_jct  are  updated,  the  corresponding  element  of  B _ is  set  to  1  and  program 

control  continues  on  to  the  simulation  stage. 

Upon  entering  the  simulation  stage,  the  matrix  B _ is  examined  to  find  which  simulations 

corresponding  to  the  range  of  control  input  magnitudes,  have  yet  to  be  run.  Then,  another 
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loop  is  started  in  which  those  simulations  are  completed.  In  this  stage^  the  following 
commands  are  executed: 

1.  Global  variables  are  initiahsed  with  ghs_i  nit 

2.  Replay  window  is  created  if  flagged  in  the  program  preferences 

3.  Trim  time-control  matrix  TOO  is  set  to  it's  initial  value 

4.  Configuration  velocity  vector  the  position  &  orientation  vector  r  and  the  full  state 
vector  X  are  set  to  their  initial  values 

5.  Time-control  input  matrix  TD  is  generated  for  the  specified  time  sequence  and  offset 
according  to  its  trim  value 

6.  Simulation  is  rim  via  hsl  SI  m 

7.  Data  struct  is  created  with  the  simulation  output^  including  time  state  controls  Q 

modified  controls  CC^  body-axes  velocities  Va^  body-axes  accelerations  Vadot^  cable 
angles  Ac load  deviation  Aj  and  cable  force  Fc 

8.  Save  the  Conf  i  g  and  Data  structs  to  file 

9.  Update  the  maximum  load  deviation  matrix  LdA 

The  last  step  of  the  simulation  stage  involves  creating  the  VRML  file  from  fhe  simulation  data 
output.  In  order  to  accomphsh  this^  firsL  the  load's  vertical  eg  location  must  be  reset  back  to 
centre^  and  the  load  position  time-history  needs  to  be  adjusted  accordingly  (in  effecL 
removing  the  offset,  Zcg),  Then,  the  user  and  fixed-inertial  viewpoint  offsets  are  defined  and 
the  VRML  file  generated  via  the  function  hs1  vrin1.  If  the  executable  GZIP  exists,  the  file  is 
compressed  from  its  native  ASCII  format  to  binary.  Lastly,  the  corresponding  element  in  the 

special  matrix  B _ is  set  to  2  to  indicate  a  successfully  simulated  flight-point  and  both  loops 

(balance/ trim  and  simulation)  are  closed.  The  progress  and/or  total  execution  time  is 
displayed  and  aU  variables  used  in  the  script  are  saved  to  file.  These  include  all  of  the  special 
variables  B_*  used  in  the  script  to  store  data  from  each  flight-point. 

If  for  any  flight-point,  the  program  has  too  much  difficulty  in  attaining  trim  (and  the  user  is 
certain  that  trim  is  possible  for  that  particular  configuration  and  flight-point),  a  different 
approach  may  be  necessary.  One  method  is  to  vary  the  initial  guess  of  the  trim  state  variables, 
rO  and  dO  to  try  and  ehminate  any  optimisation  problems  with  local  minima.  This  can  be 
done,  either  from  within  the  GUI,  or  by  executing  commands  from  the  workplace  itself.  In  the 

latter  case,  should  trim  be  attained,  it  is  imperative  that  the  special  matrix  B _ be  updated  to 

indicate  a  successfully  triormed  flight-point.  Another  alternative  is  to  modify  trim  procedure 
of  the  main  script  to  account  for  difficult  cases  or  even  just  that  particular  flight-point.  Perhaps 
the  simplest  solution  however,  is  to  extend  the  time  frame  by  some  period  preceding  the 
commanded  manoeuvre  in  which  the  heficopter  and  load  are  stabflised  using  rate- feedback. 
This  was  implemented  for  the  trim  and  simulation  of  the  CH-47D  and  Merlo  Rough  Terrain 
Forklift,  in  the  script  ch47bmrtf5l  m.  m.  The  simulation  time  was  increased  to  15  seconds,  the 
first  5  seconds  of  which,  the  heficopter  was  controlled  by  rate-feedback  on  the  stick  and 
pedals.  Consequently,  the  helicopter  is  still  able  to  stabifise  from  an  initial  state  that  is  not 
adequately  trimmed  at  several  flight-points. 
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Tliere  lure  viuious  wiiys  iii  which  the  cUita  from  a  simulation  experiment  can  be  imalysed.  Hie 
script  hslmap  is  one  example^  whicli  plots  coIollf  contour  maps  of  tlie  maximum  load 
deviation  stored  in  the  mati'ix  B_dA.  The  script  cati  be  inn  immediately  after  a  set  of 
trim /simulation  runs,  or  after  loadinf^  the  final  file  (above)  containing  all  relevant  variables 
from  such  an  experiment,  tislmap  creates  one  plot  per  control -in put  magnitude  specified  in 
the  vector  B_T  as  shown  in  Figure  3. 17^  Fig  Lire  3. 18,  and  Figure  3. 19.  hacli  plots  illustrates  the 
variation  of  maximum  load  devia  tion  with,  m  this  example,  load  mass  and  airspeed.  Contour 
lines  are  also  added  to  the  plots  for  more  clarity.  For  the  CH-47D  helicopter  and  3000 1 .  Water 
Tank  load,  it  can  be  seen  that  the  maximum  load  deviation  generally  decreases  with 
increasing  airspeed  and  load  mass,  imd  increases  with  an  increase  in  tlie  control  mput 
magnitude.  There  do  not  appear  to  be  any  anomalies  in  the  contours. 


3000L  Water  Tank  Load  (3330-9943  lb) ;  Doublet  Input  =  0.0  in  Max  Load  Dev'n  (deg) 
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rigure  3. 1 7;  Max  Load  Deuiation  Contour  Plot  CHM7D  +  3000  L  Water  Tank,  Lateral  Do^Mct 

Input  Magnitude  =  0.0  in 
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3000L  Water  T ark  Load  (3330-9943  lb) ;  Doublet  Input  =  1 .0  in  Max  Load  DeVn  (deg) 


Figure  3.18:  Max  Load  Deviation  Contour  Plot  CHA7D  +  3000  L  Water  Tank,  Lateral  Doublet 
Input  Magnitude  =1.0  in 


3000L  Water  T ark  Load  (3330-9943  lb) :  Doublet  Input  =  2.0  in  Max  Load  Dev'n  (deg) 
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Figure  3.19:  Max  Load  Demation  Con  Imr  Plot:  CH-47D  +  3000  L  Water  Tank,  Lateral  Doublet 
Input  Magnitude  =  2.0  in 
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Another  script,  hs1  rng  performs  essentially  the  same  operation,  but  presents  the  data  in  a 
different  way,  as  shown  in  Figure  3.20.  With  this  script,  the  maximum  load  deviation 
magnitude  is  plot  against  airspeed  using  lines  of  constant  load  mass.  There  is  one  group  of 
lines  for  each  control  input,  all  drawn  on  the  same  graph. 


3000L  Water  Tank  Load  (3330-9943  lb) 


Airspeed  (kt) 


Figure  3.20:  Max  Load  DeznaUon  Range  Plot:  CHA7D  +  3000L  Y^ater  Tank 


Both  hsimap  and  hsl  rng  will  work  for  most  simulation  experiments.  However,  they  may 
need  to  be  edited  if  different  data  is  to  be  examined,  or  the  resulting  data  is  not  presented  in 
the  manner  desired.  For  example,  if  the  lateral  eg  position  is  specified  as  an  independent 
variable  in  the  suite  of  flight-points,  or  if  the  user  wishes  to  fit  a  polynomial  curve  to  the  data 
prior  to  display. 

Of  course,  these  files  are  just  scripts  that  load  the  configuration  and  output  data  for  each 
flight-point  and  then  compile  the  relevant  information  into  matrices  which  are  subsequently 
displayed  in  various  formats.  For  each  experiment,  the  files  that  are  saved  include:  one 
containing  the  initial  configuration  information;  one  containing  the  configuration  information 
and  output  data  for  each  flight-point  simulated  (totalling  B_inN  *  B_VN  *  B_TN);  and  one  with 
all  of  the  special  and  extraneous  variables  used.  In  the  example  above,  these  are: 
ch47b3kltsim-config.mat,  ch47b3k1t5im-v060M333D00.mat  ...  ch47b3kltsim- 

V120M994D20 . mat  and  ch47b3k1ts'im.tnat. 

The  first  file  contains  the  Conf  i  g  struct,  the  second  through  to  the  second  last  contains  both 
Conf  ig  and  Data  structs  and  the  last  contains  all  variables  resident  following  last  iteration  of 
the  main  loop.  Therefore,  in  order  to  access  the  datafrom  any  particular  flight-point,  the  user 
needs  simply  to  load  the  file  via  the  command: 

»  load  Ch47b3k1ts'im-V090M627D20 
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which  brings  the  Conf  i  g  and  Data  stmcts  into  the  workspace.  As  detailed  below  (Appendix 
A2.1),  Data  contains  fields  time,  state,  control,  modcontro!,  bodyvel,  bodyacc, 
cableangle,  loaddeviation,  and  cableforce.  So,  to  plot  the  cable  force  variation  against 
time,  issue  the  command: 

»  plot CData. time, Data. cableforce) 

As  another  example,  the  script  ghsl_ct  rl  generates  the  time-control  input  matrix,  TD  for  a 
control  sequence  specified  in  Config,  so,  to  plot  the  control  inputs  against  time,  issue  the 
commands: 

»  N  =  config. si mdata.numsteps;  t  =  Data. time;  ghs1_ctr1 
»  for  j  s  1:4,  subp1otC2, 2, j) ,  p1otCt,TDC: , j+1,1)),  end 

Lastly,  it  is  veiy  easy  to  create  a  script  similar  to  hsl  map  and  hsl  rng  that  cycle  through  the 
flight-points,  load  the  data  from  each  one,  and  compile  specific  values  into  a  large  matrix  for 
display.  The  following  commands  do  just  that,  for  the  cable  roll  (lateral  swing)  angle  over  a 
range  of  airspeeds. 

»  B_title  =  'ch47b3k1tsim' ;  B_mdir  =  [pwd, '\HSL\' ,B_title, ; 

»  1oadC[B_mdir,B_tit1e]),  B^P  =  zerosCN,B_VN); 

»  for  B_j  =  1:b_vn 
»  B_i  =  1;  B_k  =  3; 

»  Bjnfile  =  sprintf  C[B_title, '-V%03drtlK)3dDSI502d']  ,  ... 

»  round  Cb_vCb_3))  ,  round  CB_mCB_i)/10) ,  round  Cb_tCbJs)*  10)); 

»  fprintf ([Bjnfile, '\n']),  loadC[B_mdir,B_mfil6]) 

»  B_dPC:,B_j)  =  Data.cableangleC: , 1) ; 

»  end 

»  surf CB_V,t,B_dP* 180/pi ) ;  viewC[62. 5,30]) 

»  xlabelC’ Air  speed  Ckt)'),  ylabel  Crime  (s)') 

»  title ('Cable  ro11  angle  (deg)') 

The  result  from  this  script  is  shown  in  Figure  3.21. 


Cable  rol  angle  (deg) 


AJrspeed  (kt) 


Tinte  (s) 


Figure  3.21:  Cable  Roll  Angle:  CH-47D  +  3000  L  Water  Tank,  Load  Mass  =  3330  lb;  Lateral 
Doublet  Input  Magnitude  =2.0  in 
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Full  details  of  the  data  matrices  are  outlined  in  [20].  For  a  single-cable  configuration^  the 
system  state  matrix: 


X  =  [u  if]  (4) 

where  the  configuration  velocity  and  position  matrices: 

^  =  [4  4  4  Pi  91  h  P2  92  ^2]  .c. 

r  1 

R  =  [x^  yi  zi  ^  0^  (//^  X2  72  ^2  h  ^2 

and  [j8^  ^  £[]  represent  the  hehcopter  inertial  velocities^  [jk^  3%  ^  ]  represent  the  load 

inertial  velocities^  l^i  7^ l^i  the  cable  velocities  in  cable  axes^  0i  t//i  ] 

are  the  hehcopter  Euler  angles^  and  [^2  ^2  ^2]  the  load  Euler  angles.  For  other 

configurations,  the  velocities  used  in  place  of  the  cable  velocities  are  different,  but  still  relative 
to  the  hehcopter. 

The  control  matrix: 


C  =  [Si,  S,  S,  S,]  (6) 

where  denotes  the  longitudinal  stick  control,  denotes  the  lateral  stick,  S^.  denotes  the 
pedals  and  denotes  the  collective  stick.  Also,  the  augmented  control  matrix  CC  has  the 
same  control  inputs,  modified  by  any  control-feedback  commands  present. 

The  body  velocity  matrix: 

K=W  Wi  U2  V2  W2  Pi  qi  ri  P2  q2  ?'2]  (7) 

and  the  body  acceleration  matrix  l^is  just  the  differential  of  . 

The  cable  angle  matrix: 


Aj=y>cj  ^cj  v'cji  (8) 

comprises  the  roU,  pitch  and  yaw  angles  of  the  first  cable  (forward-starboard)  supporting  the 
load.  The  load  deviation  vector  Ajj ,  represents  the  absolute  angle  made  by  a  line  from  the 

hehcopter  eg  to  the  load  eg  from  its  trim  orientation. 

Lastly,  the  cable  force,  FC  is  a  vector  of  the  non-dimensional  tensile  force  within  the  cable. 

Ah  of  fhe  data  matrices  and  vectors  wih  normally  have  the  length  specified  in  the  simulation 
configuration.  Note,  however,  that  if  the  simulation  has  diverged,  or  rim  into  error,  or 
stopped  prematurely,  the  length  of  these  will  be  less. 
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3.7.2  Extending  the  Model 

Aside  from  writing  scripts  to  rim  multiple  simiiLations  or  access  the  results  for  display^  the 
user  may  wish  to  alter  the  model  itself.  Common  tasks  might  involve  changing  the  load^s 
centre-of-gravity  location,  or  adjusting  the  Limits  imposed  during  trim.  Some  of  these  aspects 
wtil  now  be  discussed. 

In  order  to  alter  the  (j-l)th  load's  eg  position,  the  only  variables  that  need  to  be  changed  are 
those  in  the  struct  LDAT_C  j  )  -  S,  which  holds  all  the  geometric  data  for  the  load.  Specifically, 
these  are:  LDAT_Cj)  .  S- Li nks.  Data  and  LDAT_Cj) -S.  Patches.  Data  —  both  3-column 
matrices  with  (x,  y,  z)  data  for  the  links  (load  attachments)  and  the  patches  (load  surface);  their 
length  is  unimportant.  The  positions  of  the  links  and  patches  are  both  relative  to  the  load's 
centre-of-gravity,  so  the  eg  can  be  effectively  moved  by  displacing  the  links  and  patches  by 
some  amount.  For  example,  in  the  script  for  the  trim  and  simulation  of  the  CH-47D  and 
ANTPQ-36  Radar,  ch47bantpqs1  m .  m,  the  eg  is  offset  horizontally  by  an  incremental  amount, 
B_x  with  the  commands: 

LDAT_C2) -S- Links. DataC :  ,1)  =  Delinks  -  BLjc(B_i): 

LDAT-C2). S. Patches. DataC ; ,1)  =  BLpatches  -  B_jcCB^i); 

Note  that  the  links  and  patches  must  be  reset,  and  the  load  position  adjusted  to  accoimt  for  the 
offset,  before  generating  a  VRML  model,  since  this  function  assumes  that  the  eg  is  at  the  load's 
centre.  It  is  also  possible  to  change  the  location  of  load's  attachment  points  in  this  way. 

To  rearrange  the  hehcopter-load  attachment  locations,  the  cable  index  matrix  for  the 
(j-l)th  load  CDAT_(  j  )  .  i  must  be  edited.  This  is  a  2-row  matrix  with  each  column  representing 
a  different  cable  in  the  configuration.  The  top  row  Hsts  the  hehcopter  attachment  indices  and 
the  bottom  row  Hsts  the  load  attachment  indices.  The  attachment  points  are  ordered  aftward 
from  the  forward-most  location,  or  in  a  clockwise  direction  (looking  down)  from  the  forward 
starboard  location,  as  detailed  in  Section  3.3.2.  Typically, 

CDAT_Cj)-i  =  [  ci_heli_fwdLto_aft  ;  ci_loa(LcwLfrom_fwdstbd  ]; 

Thus,  for  a  four-cable  configuration,  the  attachment  points  wtil  be  ordered:  forward-starboard, 
aft-starboard,  aft-port,  forward-port.  If  those  cables  are  attached  in  a  tandem  configuration  to 
the  forward  and  aft  attachment  points  on  a  hehcopter  with  three  hooks,  the  matrix  wih  be 
equal  to  [Ip3,3,l;l,2,3r4].  If  the  aft  two  cables  are  brought  forward  to  the  middle 
hehcopter  hook,  the  matrix  wfh  become  [Ip2p2,l;lp2p3p4].  A  bifilar  configuration  wfh 
normally  use  the  same  hehcopter  hooks  (forward  and  aft)  as  a  tandem  configuration,  so  the 
matrix  would  be  numbered:  [1,3;  1,2].  Bringing  the  aft  cable  forward  would  change  this  to 
[1,2;1.2]  .  If  a  second,  multiple-point  load  were  placed  on  the  third  hook,  it  would  have 
indices  equal  to  [3 , 3 , 3 , 3 ;  1, 2  p  3 , 4] . 

Of  course,  the  properties  of  individual  cables  can  also  be  modified  by  directly  editing  any  of 
the  fields  within  the  cable  data  matrix,  CDAT_.  The  stiffness  and  damping  of  the  cables  are  held 
in  vectors  CDAT_C  j  )  .  K  and  CDAT_(  j  )  .  C,  respectively.  The  unstretched  and  current  lengths  are 
in  CDAT_Cj)  .10  and  CDAT_(j)  .1,  respectively.  However,  care  must  be  taken  if  adjusting 
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these  individually^  because  for  tandem  configurations^  the  forward  two  shngs  and  aft  two 
shng  must  have  the  same  lengths.  For  multiple-point  configurations^  all  four  stings  must  have 
the  same  length;  otherwise  the  load  cannot  be  set  to  its  initial  position^  before  balancing.  The 
cable  stretch  rates  CDAT_(  j  )  .  1  dot  and  the  cable  direction  vectors  CDAT_(  j  )  .  k_N/  should  not 
be  altered. 

The  various  fields  within  the  heticopter  and  load  data  matrices  HDAT_  and  LDAT_  can  also  be 
changed  directly^  but  would  be  wiser  to  edit  the  m-files  that  create  them.  For  the  heticopter^ 
these  files  are  h5l_ch47bdat  ■  m,  hs1_ch53ddat .  m,  and  hsl  juhlhdat .  m.  The  fields  include  the 
linearised  aerodynamic  coefficients  HDAT_,A^  heticopter  reference  area  HDAT_-R^  reference 
length  HDAT_,  as  well  as  the  mass  HDAT_.m^  inertia  HDAT_,  3^  control  timits  HDAT_,C  and 
geometric  structure  HDAT_.S.  For  the  loads,  there  is  hsl_ai  rfoil dat .m,  hs1_boxdat .m, 
h5l_cy11nderdat.m,  hs1_m11vandat.m,  hs1_p1atedat.m,  hsl.ribdat.m,  and 
hs1_sboxdat  -  m.  The  fields  created  in  each  include  the  aerodynamic  coefficients  LDAT_(  j  )  .  A, 
mass  LDAT_Cj)  -  m,  inertia  LDAT_Cj)  .  3,  geometric  structure  LDAT_(j  )  -  S.  The  only  field  that 
should  not  be  altered  is  the  load  angle  vector  LDAT_(  j  )  -  a. 

One  useful  experiment  involves  constraining  the  helicopter-load  model  in  certain  axes  during 
both  trim  and  simulation.  This  can  be  achieved  with  the  use  of  infinite  mass  and  inertial 
elements  in  the  heticopter  and  load  data  matrices.  For  example,  the  heticopter  and  load  can  be 
constrained  to  the  x-z  plane  by  setting:  niC2)  =  Inf, 

3(1,1)  =  Inf,  and  J  (3 , 3)  =  Inf  for  both  HDAT_  and  LDAT_Cj).  Further,  the  helicopter  can 
be  constrained  to  move  only  along  the  x-axis  by  also  setting:  HDAT_-mC3)  =  Inf  and 
HDAT_.  3(2,2)  =  Inf.  The  matrix  that  is  actually  used  in  the  functions  is  the  combined 
helicopter-load  mass-inertia  matrix  D,  which  is  created  in  ghsl_i  ni  t  by  concatenating  each  of 
the  sub  matrices  along  the  main  diagonal. 

Finally,  it  may  become  necessary  to  increase  the  timits  imposed  during  trim  if  the  true  angles 
associated  with  the  heticopter  and  load  orientation  are  large.  These  are  set  in  hsl  _t  ri  mf  un  -  m 
and  depend  on  the  sting  configuration  and  the  airspeed  at  trim.  The  matrix  dl  i  m  contains  the 
upper  and  lower  timits  for  the  heticopter  roll  and  pitch  angles  and  the  control  inputs.  Those 
control  timits  actually  come  from  fhe  heticopter  data  matrix  HDAT_,  C,  and  are  created  in  the 
m-files  above.  The  matrix  rl  1m  defines  the  upper  and  lower  timits  for  each  load^s  roU,  pitch 
and  yaw  angles  in  an  elastic  configuration,  while  al  i  m  defines  the  timits  for  other  load  angles 
(that  vary  according  to  the  sting  configuration)  used  in  an  inelastic  model.  Around  tiover  (< 
10ft/ s),  ttie  absolute  timits  generally  range  from  5°  to  20°,  and  above  that,  they  range  from  10° 
to  60°.  It  is  possible  to  see  the  optimisation  routine  hitting  these  timits  in  the  REPLAY  display 
during  trim  —  the  heticopter,  loads  or  controls  will  turn  red.  The  controls  are  hard-timited,  so 
if  the  timits  are  met,  ttiey  will  not  be  exceeded.  AU  of  the  angles,  on  the  other  hand,  have  only 
soft  timits  that  may  be  exceeded.  However,  any  variable  that  has  reached  or  passed  its  timit 
will  substantially  penalise  the  target  function,  which  for  trim  is  fhe  residual  (out-of-balance) 
external  force.  As  a  result,  the  optimisation  routine  should  stay  within  the  timits  specified. 
While  it  is  possible  to  maintain  realistic  limits  for  the  heticopter  and  load  angles,  it  can  be 
much  harder  to  stay  within  the  control  timits  of  the  heticopter.  If  the  controls  are  hitting  then 
timits  for  a  range  of  different  flight-points,  it  is  likely  that  the  helicopter-load  system  is  outside 
its  envelope  and  will  not  be  able  to  achieve  trim. 
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4.  Conclusion 

A  system  description  has  been  provided  for  the  hehcopter  slimg  load  simulation  software 
package  HSLSIM^  followed  by  a  detailed  User^s  Guide.  Examples  are  provided  throughout  for 
a  CH-47D  hehcopter  carrying  a  3000  L  water  tank  and  responding  to  a  lateral  doublet  control 
input.  A  description  of  the  simulation  procedure  is  detailed  for  both  GUI  and  script-driven 
execution.  The  model  should  provide  Army  personnel  with  an  insight  into  dynamic  stabhity 
of  slimg  loads  before  flight  testing  thereby  providing  improved  safety  and  reduced  flight 
testing  time. 
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Appendix  A:  MATLAB®  Toolbox 


A.l.  Package  Contents 

A.l  .1  MATLAB®  M-files,  MEX-f iles  and  M AT-files 


a2q.iii 

a1t2s1g.m 

ch47b3k1ts1m.in 

ch47bantpqs1m.m 

ch47bmrtfs1m.m 

ch47btc1s1m.m 

count,  m 

dtime.m 

e2q.m 

e2v.m 

erates.m 

eratesi .m 

euler.m 

eulers.m 

ghs1_about.m 

ghs1_ctr1 .m 

ghs1_ctr1cfg.m 

ghs1_he11cfg.m 

ghs1_1n1t.m 

ghs1_1oadcfg.m 

ghs1_p1ot.m 

ghs1_pref5.m 

ghs1_tr1m1n1t.m 

ghsisim.m 

hs1_a1 rfo1 1 dat . m 

hsljjifilar.m 

hs1_boxdat.m 

hs1_ch47bdat.iii 

hs1_ch53ddat.m 

hs1_checkopts.m 

hsl_coeffs.m 

hs1_conf1g.m 

hs1_cy1 1 nde  rdat . m 

hs1_forces.m 

hs1_1nd1ces.m 

hs1_1oadforces.m 

hs1_1oadfun.m 

hs1_in11vandat.m 

hs1_output.m 


Angular  displacement  conversion;  Angxilar^Quaternion  representation 
Altitude  to  density-ratio  conversion 
Simulation  script  for  CH-47D  +  3000L  Water-tank  system 
Simulation  script  for  CH-47D  +  ANTPQ-36  Radar  system 
Simulation  script  for  CH-47D  +  Merlo  Rough  Terrain  Forklift  system 
Simulation  script  for  CH-47D  +  Trailer  Cargo  Lightweight  system 
Graphical  progress  cormt  &  time  indicator 

Time  string  formatted  in  years,  months,  days,  hours,  mins  &  secs 
Angular  displacement  conversion;  Euler->Quaternion  representation 
Angular  displacement  conversion;  Euler->VRML  rotation  representation 
Euler  transform  matrix  for  angular  velocities 
Inverse  Euler  transform  matrix  for  angular  velocities 
Transformation  matrix  for  an  Euler-angle  triplet 
Transformation  matrix  for  an  Euler-angle  rotation  sequence 
GHSLSIM  function  to  create  'About... '  dialog  window 
GHSLSIM  function  to  generate  control  inputs  from  config  data 
GHSLSIM  function  to  create  control  input  scheme 
GHSLSIM  function  to  configure  hehcopter 

GHSLSIM  function  to  initiahse  trim/simulation  with  config  data 
GHSLSIM  function  to  configure  load(s) 

GHSLSIM  function  to  plot  output  from  simulation 
GHSLSIM  function  to  set  program  preferences 
GHSLSIM  function  to  set  initial  trim  state 

Graphical  interface  to  the  Hehcopter  Slrmg-load  Simulation  program 
Airfoil  (wing)  aerodynamic,  mass,  inertial  and  geometric  data 
Determine  load  position  for  bifilar  sling  configuration 
Box  container  aerodynamic,  mass,  inertial  and  geometric  data 
CH-47B  aircraft  aerodynamic,  mass,  inertial  and  geometric  data 
CH-53D  aircraft  aerodynamic,  mass,  inertial  and  geometric  data 
Check  input  parameters  for  HSLSIM 
Create  the  aerodynamic  state-space  matrices 

Create  the  configuration  matrices  and  the  basis  matrix  for  solution 
Cylindrical  container  aerodynamic,  mass,  inertial  and  geometric  data 
Compute  the  accelerations  for  a  hehcopter  slrmg-load  system 
Generate  subsystem  state  vector  indices  for  one  body 
Compute  the  aerodynamic  force  for  one  load 

Minimisation  function  calculates  load  position  &  orientation  for  balance 
MILVAN  container  aerodynamic,  mass,  inertial  and  geometric  data 
Output  display  function  for  HSL_FORCES 
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hs1_p1atedat.m 

hs1_r1bdat.m 

h5l_sboxdat.m 

hs1_set1oad.m 

hsl_trimfun,m 

hs1_uhlhdat.m 

hsUoad.m 

hsimap.m 

hsiplot.m 

hs1 rng.m 

hs1s1m.m 

hsitrim.m 

hs! vrm1 . m 

j2e.m 

jacoblan.m 

kZe.m 

1abe1f1x.m 
legendtrans.m 
11nspace2.m 
odef  .m 

pendfreq.m 
q2e.m 
q2v.m 
qratesi .m 
quaternion. m 
replay. m 
setportralt.m 
skews,  m 
subplot.. m 
t2e.m 
v2a.m 


Flat  plate  aerodynamic^  mass^  inertial  and  geometric  data 
Rigid  Inflatable  Boat  aerodynamic^  mass^  inertial  and  geometric  data 
Small  box  Container  aerodynamic^  mass^  inertial  and  geometric  data 
Calculate  the  load  position  from  the  hehcopter  position  and  orientation  as  well 
as  various  load  and  cable  angles 

Miniinisation  function  to  calculate  hehcopter  and  load  position  and  orientation 
for  trim 

UH-IH  aircraft  aerodynamic^  mass^  inertial  and  geometric  data 

Determine  the  load  position  through  optimisation 

Colour  map  graphs  of  max  load  deviation  from  multiple  HSLSIM  rims 

Time-history  plots  of  states  and  controls  from  simulation 

line  plots  of  maximum  load  deviation  from  multiple  HSLSIM  rims 

Hehcopter  Slung-Load  Simulation 

Determine  the  hehcopter  and  load  trim  state  through  optimisation 
Generate  dynamic  VRML  model  for  replay  of  HSLSIM  simulation 
Convert  j-unit  vectors  to  wing-oriented  Euler  angles  used  in  HSLVRML 
Numericahy  compute  the  Jacobian  dF/dY  of  function  F(T^Y) 

Convert  k-unit  vectors  to  cable-oriented  Euler  angles 
Fix  rotation  and  ahgnment  settings  for  y-labels  on  axes 
Place  transposed  legend  on  current  window 
Linearly  spaced  matrix. 

Differential  equation  solver  based  on  Runge-Kutta  integration  with  fixed  step- 
size 

Estimate  slung  load  frequencies  using  equivalent  pendulum  system 
Angular  displacement  conversion;  Quaternion  ->  Euler  representation 
Angular  displacement  conversion;  Quaternion  ->  VRML  rotation 
Inverse  quaternion  transform  matrix  for  angular  velocities 
Transformation  matrix  for  a  quaternion  set 
Display  dynamic  system  response  path 
Set  current  paper  positions  for  portrait  orientation 
General  skew-symmetric  matrix  S(x^y^z) 

Create  axes  in  tiled  positions  with  strict  scaling 
Convert  transformation  matrix  to  Euler  angles 
Calculate  quaternion  transform  defined  by  two  vectors 
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mexrotorgen.d11 


Compute  aerodynamic  forces  using  RotorGen  MEX  model 


ghs1 default. mat 
ghs1demo_*.mat 
ghs1 demo_3  k1 1 . mat 
ghs1  demo_3  k1  ^t  r  1  m .  mat 
ghs1  demo_3k1  t^run .  mat 
ghs1 demo_b1 f 1 1 ar . mat 
ghs1 demo_mu1 t1 load . mat 
ghs1 demo_mu1 ti p1 e . mat 
ghs1 demo_no1 oad . mat 
ghs1 demo_si ng1 e . mat 
ghs1 demo_s1 ng1 e_t r1 m . mat 
ghs1 demo_s1 ng1 e_run . mat 
ghs1 demo_tandem. mat 
ghs1 demo_tandem_t r1 m . mat 
ghs1 demo_tandem_run . mat 
ghsiprefs.mat 


Default  configuration  file 

Demonstration  configuration  flies  with  the  CH-47D  heflcopter 
3000L  Water-tank  load;  Multiple-point  sflngs 
As  above  with  stable  trim  estimate 
As  above  with  simulation  data 
Box  container;  Bifllar  sflngs 

Two  Box  containers;  Single  and  multiple- point  sflngs 
Box  container;  Multiple- point  sflngs 
No  load  (heflcopter  only) 

Box  container;  Single-point  sling 
As  above  with  stable  trim  estimate 
As  above  with  simulation  data 
Box  container;  Tandem  sflngs 
As  above  with  stable  trim  estimate 
As  above  with  simulation  data 
Program  preferences  file 


A.1.2  VRML  Model  Files 


antpq.wr1 

ANTPQ-36  Radar  load 

box.wr1 

Generic  rectangular  load 

ch47b.wr1 

CH-47B  Chinook  Hehcopter  with  two  main  rotors 

ch47b-tex.wri 

As  above  with  image-based  texture  mapping 

cyl1nder.wri 

Generic  cylindrical  load 

ground. wr1 

Groimd  terrain  mesh  with  random  height  and  texture  mapping 

mms.wr1 

MMS  Medium  Maintenance  Shelter  load 

mms-tex.wr1 

As  above  with  image-based  texture  mapping 

p1ate.wr1 

Generic  plate  load 

rib.wrl 

RIB  Rigid  Inflatable  Boat  load 

s70a9.wr1 

S70- A9  Blackhawk  Hehcopter  with  main  and  tad  rotors 

tcl  .wrl 

TCL  Trailer  Cargo  Lightweight  load 

uhlh.wr1 

UH-IH  Iroqois  Hehcopter  with  main  and  tail  rotors 
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A.2.  Code  Listing:  Multiple  Simulations 


A. 2.1  ch47b3kltsim.m 


CH47B3KLTSIM  :  Simulation  script  for  CH-47D  +  3000L  Water-tank  system 
CH47B3KLTSIM 
Edi t  for  detail s  . . . 

R.A.  Stuckey  20/01/03  (c)  2003,  Defence  Science  and  Technology  Organisation 


%  CONFIGURATION  SECTION  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  Check  if  script  has  been  run  before  during  the  same  session 

if  -existC'B _ ','var') 

%  First,  the  title  and  directory  used  for  all  saved  files  must  be  set 
B_tit1e  =  'ch47b3k1tsim' ; 

B_mdir  =  [pwd , ' \hsl\ ' , B_ti t1 e , ' \ ; 
if  ~exi St CB_mdi r) ,  mkdi r CB_mdi r) ;  end 

%  As  well  as  the  range  of  load  masses  used 
B_m  =  1  inspace  (3330, 9943 , 10) : 

%  And  the  range  of  airspeeds  at  which  the  helicopter  is  trimmed 
B_V  =  [60:10:120] ; 

%  Also,  the  control  input  magnitudes 
B_T  =  [0.0,1.0,2.01; 

%  Next,  the  preferences  associated  with  ghsisim  need  to  be  specified 
Pref s  =  struct C  . . . 

' webbrowser ' , ' ' ,  ... 

'vrmidir',  [pwd,  ' \- ■  \VRML\HSLy  1 ,  ■■■ 

' gzi pexe c:\Program  Fi1es\winzip\gzip.exe',  ... 

' replaywi n ' , 1,  ... 

' vtrackwi n ' , 1,  ... 

' countwin ' , 1,  ... 

' verbose ' , 1  ... 

): 


Create  the  configuration  struct,  which  contains  all  of  the  remaining  data  for 
the  simulation,  helicopter,  load,  controls  and  trim. 

Config  =  structC 


'  simdata' , structC 
' intti me ' , 

' numsteps ' , 

' si  maxes ' , 

'  angrep'  , 

' si mtype ' , 

' si ingmod ' , 

' trimf un ' , 

' si mfun ' , 

' trimvel ' , 

'  trimalt '  , 

' numloads ' , 


0.100,  ... 

101,  ... 

'  combi ned ' , 

' euler ' ,  .  . 

'  noni i near' 

' el astic '  , 

'  1 sqnoni in' 

'  ode45 ' ,  . . 
[B_VCl)/0. 5925,0,0], 
100.0,  ... 

1  .  .  . 


j ,  ... 

' hel idata ', StructC  ■■■ 

' hel itype ' ,  ' ch47' ,  ... 

' si mmodel ' ,  ' mex'  ,  ... 

' seas ' ,  0,  .  .  . 

' hel iwei ght ' ,  46000,  ... 

' hel i inertia' ,  [43200, 0, 0; 0, 250000, 0; 0, 0, 234000]  . . . 

).  ... 

' loaddata ', StructC  ... 

' loadtype ' ,  { ' box' },  ... 

' loadwei ght ' ,  {B_mCl)},  ... 

' loadi nerti a' ,  {diagC[  0.33  1.20  1.20  ]*B_mCl))},  ... 

'loadsize',  {[96.0,77.5,46.01/12},  ... 

'slingefg',  {'multiple'},  ... 

' cablelength ' ,  {16.0  +  3.5},  ... 

' cablestiff ' ,  {20000.0},  ... 

'cabledamp',  {500.0},  ... 

'attachoffset' , { [1, 1, -1; -1, 1, -1; -1, -1, -1; 1, -1, -l]*diagC[96 . 0, 77.5,46.01/2/12)}, 
' cables] ack ' ,  {0}  . . . 

).  ... 

'  ctrldata ' , structC  . . . 
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): 


'  inputtype ' , 

'  startti me ' , 

' ramptime ' , 

' holdtime ' , 

' magni tude ' , 

' rategai n ' , 

),  ... 

' tri mdata ' , structC 
'  rO'  , 

'dO'  , 

'  xt '  , 

)  ... 


{{ ' ratefb ' } , { ' doublet ' , ' pul se ' } , { ' doublet ' } , { ' none ' }} ,  . . . 
{{0},{1.0,4.5},{1.0},{0}},  ... 

{{0},{0.5,0.5},{0.5},{0}},  ... 

{{10}, {1.0, 0.5}, {1.0}, {0}},  ... 
{{0},{B_TC1),-B_TC1)},{B_TC1)^^0.5},{0}},  ... 
{{40},{0,0},{0},{0}}  ... 


[6^0, 0,0, 0,0,0, 5*pi/lS0,0,0,0,0] '  ,  .  .  . 
[0,0,0,4.56],  ... 
zeros (1, 12)  .  .  . 


save  C [B_mdi r , B_ti tl e , ' -conf i g ' ] , ' Conf i g ' ) 

%  Lastly,  specify  the  VRML  geometric  helicopter  and  load  models  to  be  utilised 
wrlfiles  =  {  [],  [Prefs.  vrmldi  r,  '  ch47b.  wrl  '  ] ,  [Prefs .  vrml  di  r,  '  box.  wrl  '  ]  }; 

%  END  CONFIGURATION  SECTION  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

global  opt_  pref_ 

%  Create  the  main  index  matrix  which  specifies  the  order  of  simulation  runs 
B_mN  =  length CB_m);  B_VN  =  length Cb_v);  b_tn  =  length Cb_t); 

reshape C[l:B_mN*2'''cei  1  Cb_VN/2)]  ,  B_mN'''2,  cei  1  Cb_VN/2))  ;  ans  C[l:B_mN,  B_mN'^'2  :  -1:  B_mN+l] ,  : ) ; 
B_l  =  reshape  Cans  Cl  :B_mN^^B_VN),B_mN,B_VN);  B_l2  =  zerosCB_mN,  B_VN) ; 

%  Create  some  more  special  variables  used  in  the  main  loop  of  the  script 
n  =  confi g. simdata. numloads  +  1; 

B_dA  =  zeros CB_mN,B_VN,B_TN);  B _  =  B_dA; 

B_u0  =  zerosCn*6 , B_mN, B_VN) ;  B_r0  =  B_u0;  B_xt  =  B_u0;  B_d0  =  zerosC4 , B_mN, B_VN) ; 
[B_i,B_j]  =  findCB_i  ==  1); 

%  Set  the  position  &  orientation  vector,  the  control  vector  and  the  trim  state 
%  vector  for  the  first  step 

B_r0C: ,  B_i ,  B_j  )  =  Config  .  trimdata.  rO; 

B_d0C: , B_i , B_j )  =  Config . trimdata. dO' ; 

B_xt C: , B_i , B_j )  =  Config . trimdata. xt' ; 

el  se 

fprintfC'\n  Resuming  script  ...\n') 


%  Set  some  global  s  used  by  hsltrim  and  hslsim  in  each  simulation  run 
global  HDAT_  LDAT_  CDAT_ 

B_TR  =  cputime; 

%  initialise  the  hsltrim  variables 
ghsl_ini  t 

B_links  =  ldat_C2) . s. Links. DataC : , 3) ; 

B_patches  =  LDAT_C2) . s. Patches . DataC: , 3) ; 

%  start  the  main  loop 
for  B_ii  =  l:B_mN*B_VN 

%  Identify  the  current  airspeed  and  load  mass  vector  indices 
[B_i,B_j]  =  findCCB_l  ==  B_ii)&CB _ C <  2)); 

%  Skip  this  step  if  the  model  has  not  been  trimmed  and  simulated  at  this  point 
if  ~i semptyCB_i ) 

%  Adjust  the  load  mass  and  inertia,  and  the  trim  velocity  for  this  step 
Conf ig . loaddataCl) . loadwei ght  =  B_mCB_i); 

conf ig . loaddataCl) . loadi nerti a  =  diagCE  0.33  1.20  1.20  ]*B_mCB_i)); 

Conf ig . si mdata. tri mvel Cl)  =  B_vCB_j)/0. 5925; 

%  Display  the  current  trim  step 

if  Cb _ CB_i  ,B_j,  1)  <  1) 

f printf C ' v  =  %5.1fkn  m  =  %4.0flb 

' ,  B_vCB_j)  ,  B_mCB_i  ))  ; 

end 

%  Set  the  initial  trim  index  either  to  the  current,  or  to  the  previous  value 
i  f  B_i  i  ==  1 

B_ii2  =  B_ii;  B_i2  =  B_i ;  B_j2  =  B_j  ;  B_l2CB_i , B_j)  =  B_ii2: 
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el  se 

B_ii2  =  [B_i2,B_j2]  =  findCB_l  ==  B_ii2);  B_l2CB_i  , B_j)  =  B_ii2; 

end 

%  Loop  while  the  system  does  not  trim 
while  Cb _ CB_i,B_j,l)  <  1) 

%  Each  time,  setting  the  initial  trim  variables  to  previous  values 
config.trimdata.ro  =  B_rOC : , B_i 2, B_j 2) ; 
config .  trimdata.  do  =  B_dOC : ,  B_i  2,  B_j  2y  ; 
config . trimdata. xt  =  B_xtC : , B_i 2, B_j 2) ' ; 

%  BALANCE  loads:  INELASTIC  CONFIGURATION  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  Set  the  sling  model  to  inelastic 

Config . simdata. si i ngmod  =  'inelastic'; 

%  Initialise  the  hsltrim  variables 
ghsl_i  ni  t 

%  Calculate  and  adjust  the  load's  vertical  eg  location 

zcg  =  (25.0  -  C23.0*B_mCl)  +  23 . 0*CB_mCB_j)-B_mCl))A2/(B_mCend)-B_mCl)))/B_mCB_j))/12 ; 
LDAT_C2) . S. Li nks. DataC: , 3)  =  B_links  -  zcg; 

LDAT_C2) . S. Patches . DataC: , 3)  =  B_patches  -  zcg; 

if  pref_. verbose,  fprintfC'\n  Balancing  loads  (inelastic  config)  ...  \n'),  end 

%  set  the  initial  load  position  directly  beneath  the  fixed  helicopter 
jj  =  3+[l:3]; 
if  ~any(r0([jj,n*3+jj])) 

f0( [j  j  ,  n*3+j  j])  =  hsl_setload(uO,  rO,  [] ,  j  ) ; 

end 

%  Bring  up  the  replay  window  if  requested 
i  f  Prefs .  repl aywin 

hr  =  replay(0,  rO(l  :n*3)  '  -r0(  [1 :  3]  '  *ones(l,  n)  ,  1)  '  ,  r0(n*3+[l  :n*3]  )  '  ,  []  ,  S)  ; 

end 

%  Balance  the  load  beneath  the  fixed  helicopter 
hsl  load 

%  Update  the  load  positional  trim  data 
Config.trimdata.ro  =  rO; 

%  BALANCE  loads:  ELASTIC  CONFIGURATION  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  Set  the  sling  model  to  elastic 

Config . simdata. si i ngmod  =  'elastic'; 

%  Initialise  the  hsltrim  variables 
ghsl_i  ni  t 

if  pref_. verbose,  fprintf('\n  Balancing  loads  (elastic  config)  ...  \n'),  end 

%  Set  the  initial  load  position  directly  beneath  the  fixed  helicopter 
jj  =  3+[l:3]; 
if  ~any(rO([j  j  ,  n*3+jj])) 

r’0(  [j  j  ,  n*3+j  jl)  =  hsl_setload(uO,  rO,  []  ,  j  ) ; 

end 

%  Bring  up  the  replay  window  if  requested 
i  f  Prefs .  repl  aywin 

hr  =  replay(0,  rO(l  :n*3)  '  -r0(  [1 :  3]  '  *ones(l,  n) ,  1)  '  ,  r0(n*3+[l  :n*3] )  '  ,  [] ,  S) ; 

end 

%  Balance  the  load  beneath  the  fixed  helicopter 
hsl  load 

%  update  the  load  positional  trim  data 
Config.trimdata.ro  =  rO; 

%  TRIM  HELICOPTER  +  LOADS  (ELASTIC  CONFIG)  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  initialise  the  hsltrim  variables 
ghsl_i  ni  t 

%  Bring  up  the  replay  window  if  requested 
i  f  Prefs .  repl  aywin 

hr  =  replay(0,  rO(l  :n*3)  '  -rO(  [1 :  3]  '  *ones(l,  n)  ,  1)  '  ,  rO(n*3+[l  :n*3]  )  '  ,  []  ,  S)  ; 

end 
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Set  the  trim  time-control  matrix  to  it's  initial  value 
TDO  =  zerosCl,  5,2);  100(1,2:5,1)  =  dO; 

Trim  the  helicopter  load  system 
hsi  tri  m 


update  the  trim  control  vector 
do  =  TDOCl,  2:  5, 1); 

Check  the  residual  (out-of-balance)  external  force  sum  against  some  tolerance 
if  (ff  >  le-1) 

Trim  point  not  reached 
if  B_ii2  ==  1 

No  earlier  trim  points  remaining;  Exit  the  script 
fprintf  ('  \n  '  ) 

error  (sprintfC  Large  error:  F  =  %g  N  =  %g  ;  Exi  ting  .  .  .  '  ,  ff ,  Bi  i ) )  ; 
el  se 

Try  earlier  trim  point  for  initial  values 

B_ii2  =  B_ii2-1;  [B_i2,B_j2]  =  find(B_l  ==  B_ii2);  B_l2(B_i , B_j )  =  B_ii2; 


fprintf  ('  \n '  ) 

warni ng (spri ntfC  Large  error:  F  =  %g  N  = 

%4 .  Ofib  '  ,  ff ,  B_i  i ,  B_V(B_j  2) ,  B_m(B_i  2))) ; 
fprintf  ('  \n '  ) 

end 

%  Flag  the  corresponding  special  matrix  element  to  denote  an  untrimmed  state 
B_(B_i  ,B_j,  :)  =  0; 
el  se 

%  Trim  point  successfully  reached  -  update  the  special  matrices  for  the 
%  configuration  velocity,  position  &  orientation,  controls  and  trim  state 
B_uO( : , B_i , B_j)  =  uO;  B_rO(: , B_i , B_j)  =  rO; 

B_dO(: ,B_i ,B_j)  =  dO';  B_xt( : , B_i , B_j)  =  XT'; 

config.trimdata.ro  =  rO; 

Confi g. tri mdata. do  =  dO; 

Confi g. tri mdata. xt  =  XT; 

%  Flag  the  corresponding  special  matrix  element  to  denote  a  trimmed  state 
B _ (B_i , B_j , :)  =  1; 

end 

end 

%  Find  the  simulation  cases  that  have  yet  to  be  run  for  the  current  trim  point 
B_K  =  find(B _ (B_i,B_j,:)  <  2)'; 

%  And  loop  through  those 
for  B_k  =  B_K 

%  update  the  control  magnitudes 

config . Ctrl  data (2) . magnitude  =  {B_T(B_k) , -B_T(B_k)} ; 
config . Ctrl data(3) . magnitude  =  {B_T(B_k)*0. 5}; 

%  Display  the  current  simulation  step 

f  printf  ( '  D  =  %3  .  If  i  n\n '  ,  B_T(B_k))  ; 

%  SIMULATE  HELICOPTER  +  LOADS  (ELASTIC  CONFIG) 

%  Initialise  the  hsisim  variables 
ghs1_i  ni  t 

%  Bring  up  the  replay  window  if  requested 
i  f  Prefs .  repi aywin 

hr  =  replay(0,  rO(l  :n*3)  ' -r0(  [1 :  3]  ' *ones(l,  n) ,  1)  ' 

end 

%  Set  the  trim  time-control  matrix  to  its  initial  value 
TDO  =  zerosd,  5,2);  TD0(1, 2:5,1)  =  dO; 

%  Set  the  configuration  velocity  vector,  the  position  &  orientation  vector  and 
%  the  full  state  vector  to  their  initial  values 
u  =  uO;  r  =  rO;  x  =  [u;  r] ; 

%  Generate  the  control  inputs  for  the  specified  time  sequence 
ghsl_ctrl 

%  offset  the  time-control  matrix  by  its  trim  value 

TD(:, 2:5,1)  =  TD( :  ,  2 :  5 ,  l)+ones  (N ,  1)^^TD0(1 , 2  :  5 , 1)  ; 

%  Run  the  simulation 


;  Attempting  re-trim  from  V  =  %5.1fkn  m  = 


,  rO(n*3+[l:n*3])' , [],S); 
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hs1  sim 

%  dose  the  replay  window  if  present 
i  f  Prefs .  repi aywin 
close  (hr) 

end 

%  Extract  the  (open  and  closed  loop)  time-control  input  matrices 
c  =  tdC:, 2:5,1);  cc  =  tddC:,2:5); 

%  create  a  struct  with  the  simulation  output 
Data  =  struct C time t,  ... 

' state ' ,x,  ... 

' control ' , C,  ... 

' modcontrol ' , CC,  . . . 

' bodyvel ' , Va,  ... 

' bodyacc ' , Vadot,  . . . 

' cableangle ' , Acj ,  ... 

' loaddevi ati on ' , Aj j ,  ... 

' cableforce ' , FC) ; 

%  Append  the  m-file  name  with  airspeed,  load  mass  and  control  magnitude  fields 
B_mfile  =  sprintf  C [B_ti  tie ,  ' - 

V%03dM%03dD%02d ' ] , round Cb_V CB_j)) , roundCB_mCB_i )/10) , roundCB_TCB_k)*10)) ; 

%  Save  the  Config  and  Data  structs  to  file 

saveC[B_mdi r, B_mfi le] , ' confi  g ' , ' Data ' ) 

%  Update  the  maximum  load  deviation  special  matrix 
B_dACB_i  ,  B_j  ,  B_k)  =  maxCabsCAj  j  ,  2)))  ; 

%  GENERATE  VRML  MODEL  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  Reset  the  load's  vertical  eg  location  to  centre  -  necessary  for  hslvrml 
ldat_C2) . s. Li nks. DataC: , 3)  =  B_links; 
ldat_C2) . s. Patches . DataC: , 3)  =  B_patches; 

%  Adjust  the  load  position  history  to  account  for  the  new  eg  location 
for  i  =  1:N 

xCi  ,  n*6+3+  [1:3])  =  xCi ,  n*6+3+[l:  3])-  [0, 0,  zcg]*eulerCxCi  ,  n*9+3+[l:  3] )) ; 

end 

%  Fill  the  first  element  of  wrlfiles  with  the  output  filename 

wrlfiles{l}  =  [Prefs.  vrmldi  r, B_ti  tie ,  B_mfi le ,'.  wrl '] ; 

%  Define  the  user  and  fixed-inertial  viewpoint  offsets 
rv  =  [  75  75  0  ]; 

dx  =  X  -  ones (size Cx, 1) , l)*xCl, :); 

%  Place  fixed-inertial  at  midpoint  +  offset,  but  above  ground 

ri  =  (mi nCdxC : , n*6+[l : 3]))  +  maxCdxC : , n*6+[l: 3])))/2  +  [  0  100  100  ]; 
if  riC3)  >  C-xCl,n*6+3)  -  5),  ri  (3)  =  -xCl,n*6+3)  -  5;  end 

if  pref_. verbose,  fprintfC  Generating  VRML  ...  \n\n'),  end 

%  Generate  the  VRML  model 

hsl  vrml  (wrl  files{:},t,xC:,  n*6+[l:  n*6])  ,  c,  1,  pi/4,  rv,  ri  ) 

%  If  gzip  exists,  compress  the  VRML  file 
if  ~i sempty (Prefs. gzi pexe) 

[s,w]  =  systemC ['"', Prefs . gzi pexe , . , Prefs . vrml di r, B_ti tl e, ' \' , B_mfil e, ' .wrl "  -c  > 

" ' , Prefs. vrmldi  r, B_ti tie , ' \ ' , B_mfil e, ' . wrz ' ]) ; 

if  isemptyCw),  del eteC[Prefs . vrml di r, B_titl e, ' \' , B_mfil e, ' .wrl ']) ;  end 

end 

%  Adjust  the  load's  vertical  eg  location  again 

LDAT_C2) . S. Li nks. DataC: , 3)  =  B_links  -  zcg; 

LDAT_C2) . S. Patches . DataC: , 3)  =  B_patches  -  zcg; 

%  Flag  the  corresponding  special  matrix  element  to  denote  a  simulated  state 
B _ CB_i , B_j ,B_k)  =  2; 

end 

%  Append  the  eputime 

B_TR  =  [B_TR,  eputi  me] ; 

if  B_ii  <  B_mN*B_VN 

%  Estimate  and  display  the  time  remaining  until  completion 

fprintfC  Progress:  %g/%g  Time  remaining:  %s\n\n',  ... 

B_i i , B_mN*B_VN , dtimeCCB_mN*B_VN-B_i i )*mean (diff Cb_tr)))) 

el  se 
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%  Display  the  total  execution  time 

fprintfC  Total  execution  time:  %s\n\n  ' ,  dti  meCB_TRCend)-B_TRCl))) 

end 

end 

end 

%  Save  all  variables,  including  the  special  variables  to  file 
save  C  [B_md  i  r ,  B_ti  t1  e  ] ) 


49 


Page  classification:  UNCLASSIFIED 


DEFENCE  SCIENCE  AND  TECHNOLOGY  ORGANISATION 
DOCUMENT  CONTROL  DATA 


1.  PRIVACY  MARKING/ CAVEAT  (OF  DOCUMENT) 


2.  TITLE 

Helicopter  SIung-Load  Simulation  Toolbox  for  use  with  MATLAB® 


3.  SECURITY  CLASSIFICATION  (FOR  UNCLASSIFIED  REPORTS 
THAT  ARE  LIMITED  RELEASE  USE  (L)  NEXT  TO  DOCUMENT 
CLASSIFICATION) 


Document 

Title 

Abstract 


(U) 

(U) 

(U) 


4.  AUTHOR(S) 
Roger  A.  Stuckey 


5.  CORPORATE  AUTHOR 

DSTO  Defence  Science  and  Technology  Organisation 
506  Lorimer  St 

Fishermans  Bend  Victoria  3207  Australia 


6a.  DSTO  NUMBER 
DSTO-TN-0855 


6b.  AR  NUMBER 
AR-014-325 


6c.  TYPE  OF  REPORT 
Technical  Note 


7.  DOCUMENT  DATE 
August  2008 


8.  FILE  NUMBER 

9.  TASK  NUMBER 

10.  TASK  SPONSOR 

11.  NO.  OF  PAGES 

12.  NO.  OF  REFERENCES 

2003/70263 

ARM  07/038 

ARMY 

49 

20 

13.  URL  on  the  World  Wide  Web 

http:  /  /  WWW.  d  sto.  d  efence.  gov.  au  /  corporate  /reports/  DST  O- 
TN-0855.pdf 


14.  RELEASE  AUTHORITY 
Chief,  Air  Operations  Division 


15.  SECONDARY  RELEASE  STATEMENT  OF  THIS  DOCUMENT 

Approved  for  public  release 


OVERSEAS  ENQUIRIES  OUTSIDE  STATED  LIMITATIONS  SHOULD  BE  REFERRED  THROUGH  DOCLJMENT  EXCHANGE,  PO  BOX  1500,  EDINBURGH,  SA  5111 

16.  DELIBERATE  ANNOUNCEMENT 


No  Limitations 


17.  CITATION  IN  OTHER  DOCUMENTS _ Yes _ 

18.  DSTO  RESEARCH  LIBRARY  THESAURUS  http:/  / web-vic.dsto.defence.gov.au/ workareas/ library  / resources/ dsto  thesaurus. shtml 

Chinook  helicopter;  Slung  load;  Modelling;  Simulation 

19.  ABSTRACT 

This  document  outlines  the  Helicopter  Slung  Load  Simulation  (HSLSIM)  Toolbox  which  is  a  set  of  utilities  for  the  simulation,  analysis 
and  display  of  the  flight-dynamic  response  of  helicopters  with  various  external  load  configurations  within  the  MATLAB®  software 
environment.  Instructions  and  examples  are  provided  for  its  operation  and  subsequent  modifications.  The  helicopter  studied  is  the  CH- 
47D,  with  the  load  types  including  rectangular  and  cylindrical  containers  as  well  as  plate  and  airfoil  shapes. _ 

Page  classification:  UNCLASSIFIED 


